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This  report  was  prepared  by  the  Metals  Behavior  Branch,  Metals 
and  Ceramics  Division,  Materials  Laboratory,  Air  Force  Wright  Aero¬ 
nautical  Laboratories,  Wright-Patterson  Air  Force  Base,  Ohio.  The 
research  reported  herein  was  conducted  under  Project  No.  2307,  "Solid 
Mechanics,"  Task  2307P102,  "Failure  Prediction  in  Metals."  Professor 
J.  N.  Yang  was  engaged  as  a  visiting  scientist  from  the  School  of 
Engineering  and  Applied  Science,  The  George  Washington  University, 
VJashington,  D.C.,  20052. 

The  idea  of  multiple  inspection  procedures  and  its  potential 
payoffs  was  proposed  by  Dr.  W.  H.  Reimann.  Valuable  and  helpful  dis¬ 
cussions  with  Messrs.  T.  Cooper  and  D.  M.  Forney  and  with  Drs.  T.  Nicholas 
and  J.  Moyzis  are  gratefully  acknowledged. 

This  report  covers  research  conducted  during  the  period  February 
1982  to  August  1982. 
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SECTION  I 
INTRODUCTION 

Current  nondestructive  evaluation  techniques  are  not  capable  of 
repeatedly  producing  correct  indications  when  applied  to  flaws  of  the 
same  length.  The  chance  of  detecting  a  given  crack  length  depends  on 
many  factors,  such  as  the  location,  orientation  and  shape  of  the  flaw, 
materials,  inspectors,  inspection  environments,  etc.  As  a  result,  the 
probability  of  detection  (POD)  for  all  cracks  of  a  given  length  has  been 
used  in  the  literature  to  define  the  capability  of  a  particular  NDE  system 
in  a  given  environment.  Some  POD  curves  are  shown  in  Figure  1  for  various 
laboratory  inspection  techniques.  Many  other  POD  curves  can  be  found, 
for  instance,  in  References  1  through  3. 

In  practical  applications,  a  nondestructive  inspection  limit, 
is  usually  specified,  which  is  a  crack  length  corresponding  to  a  high 
detection  probability  and  a  high  confidence  level.  For  instance,  the 
damage  tolerant  specification  (References  4  and  5)  requires  that  the  NDE 
system  should  be  capable  of  detecting  the  specified  crack  length,  a^^p^, 
with  a  ninety  per  cent  detection  probability  and  a  ninety-five  per  cent 
confidence  level.  The  fracture  mechanics  residual  life,  N^,  is  the  life 
for  the  crack  length,  a^p^,  to  propagate  to  the  critical  crack  length, 
a^.,  under  expected  usage  environments,  as  shown  in  Figure  2(a).  The 
return  to  service  interval,  denoted  by  N^,  is  equal  to  N^  divided  by  a 
safety  factor,  S^,  i.e.,  Np  =  N^/Sf.  If  no  crack  is  detected  during 
inspection,  the  component  is  returned  to  service  and  the  crack  length  in 
that  component  is  assumed  to  be  equal  to  a^p^,  as  shown  in  Figure  2(b). 
Then,  in  the  analysis  of  that  component,  the  crack  length  is  reset  to  be 
equal  to  a^^p^.  Therefore,  the  inspection  limit,  aj^p^,  is  also  referred 
to  as  the  reset  crack  length.  In  the  damage  tolerant  analysis,  a  safety 
factor  of  2.0  has  been  used. 

It  follows  from  Figure  1  that  there  are  two  possible  errors  that 
can  occur  in  any  inspection  situation:  failure  to  give  a  positive 
indication  in  the  presence  of  a  crack  whose  length  is  greater  than  a^^p^ 
(Type  I  error),  and  the  giving  of  a  positive  indication  when  the  crack 
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length  is  smaller  than  a^^Q^  (Type  II  error).  For  safety  critical  compo¬ 
nents  in  airframe  structures.  Type  I  error  is  of  primary  concern.  In 
the  Retirement-For-Cause  (RFC)  analysis  of  gas  turbine  engine  components, 
however,  both  Type  I  and  Type  II  errors  are  important,  because  the 
criterion  used  in  RFC  analyses  is  the  minimization  of  the  life  cycle 
cost  (LCC)  for  engine  components  (References  6  and  7).  As  a  result,  the 
reduction  for  both  types  of  error  is  one  of  the  main  objectives  of  the 
present  study. 

The  Type  I  error  allows  the  components  containing  a  crack  length 
longer  than  a^^Q^  to  return  to  service,  thus  greatly  increasing  the 
potential  safety  hazard.  For  a  given  NDE  system,  the  Type  I  error  can 
be  made  as  small  as  possible,  by  choosing  a  large  value  for  the  inspection 
limit,  However,  as  the  value  of  aj^p^  increases,  the  return  to 

service  interval,  Np,  reduces,  thus  increasing  the  frequency  or  cost  of 
inspections.  Moreover,  while  the  Type  I  error  can  be  reduced  by  increas¬ 
ing  the  value  of  the  Type  II  error  increases  accordingly,  as  shown 

in  Figures  3(a)  and  3(b). 

The  Type  II  error  rejects  good  components  and  therefore  has  an 
adverse  effect  on  the  cost  of  replacement  and  the  life  cycle  cost.  For 
a  given  NDE  system  with  a  single  inspection,  it  is  impossible  to  reduce 
the  Type  II  error  without  increasing  the  Type  I  error,  and  vice  versa. 

It  is  obvious,  therefore,  that  the  ideal  inspection  capability  of 
an  NOE  system  is  a  unit  step  function,  as  shown  in  Figure  3(c),  in  which 
both  Type  I  and  Type  II  errors  are  zero.  Unfortunately,  such  an  ideal 
NDE  system  may  be  far  from  reality. 

There  are  situations  where  the  required  a^^Q^  is  so  small  the  NDE 
system  may  not  be  able  to  detect  it  with  high  levels  of  reliability  and 
confidence.  Such  a  difficulty  may  be  circumvented  by  the  use  of  multiple 
inspection  procedures.  The  possible  improvements  for  the  inspection 
reliability  of  an  NOE  system,  through  multiple  inspection  procedures, 
will  be  demonstrated  in  this  report. 

The  POD  curve  of  an  NOE  system  is  dependent  upon  the  specific 
location  and  material  being  inspected.  Because  an  engine  disk  contains 
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a  large  number  of  critical  locations,  each  with  Its  own  POO  curve,  such 
as  bolt  holes,  cooling  holes,  rim  holes,  etc.,  the  Inspection  reliability 
for  a  disk  differs  significantly  from  that  for  a  hole,  since  a  disk  Is 
retired  If  at  least  one  hole  Is  rejected  by  the  NOE  system.  An  explor¬ 
atory  study  for  the  Inspection  reliability  of  engine  disks  and  the  effect 
of  multiple  Inspections  Is  made  herein. 

The  objectives  of  this  study  are:  (1)  to  formulate  and  derive 
mathematically  the  resulting  POD  curve  for  components  under  multiple 
Inspections,  (11)  to  investigate  quantitatively  the  potential  benefit  of 
multiple  Inspection  procedures,  (111)  to  establish  the  direction  In  which 
the  capability  of  an  NOE  system  should  be  Improved,  (1v)  to  establish 

A 

the  strategy  and  sequence  for  multiple  Inspections  to  reduce  either  Type 
I  error  or  Type  II  error,  or  both,  (v)  to  determine  the  POD  curve  for  a 
disk  from  the  POD  of  an  NDE  system,  (v1)  to  determine  the  Inspection 
reliability  of  engine  disks,  and  (v11)  to  evaluate  the  benefit  of  multiple 
Inspections  for  engine  disks.  Numerical  results  are  presented  to  demon¬ 
strate  the  potential  benefit  of  the  proposed  multiple  inspection  procedures 
for  gas  turbine  engine  disks. 
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SECTION  II 

THEORY  OF  MULTIPLE  INSPECTIONS 

In  the  literature,  the  POD  curve  for  a  particular  NOE  system  has 
been  established  from  repeated  inspection  data  in  which  the  inspector 
may  or  may  not  be  the  same.  Then  the  inspection  data  set  is  analyzed 
using  either  the  binomial  method  (e.g..  References  1-3)  or  the  method  of 
regression  analysis  (Reference  8).  It  should  be  emphasized,  however, 
that  the  binomial  and  regression  methods  assume  explicitly  that  the  result 
of  each  inspection  using  the  same  NOE  system  (whether  the  same  inspector 
or  not)  is  statistically  independent  of  the  other  (i.e.,  independent 
sampling).  Thus  all  the  POO  curves  available  in  the  literature  are  based 
on  the  premise  that  the  results  of  multiple  inspections  are  statistically 
independent.  Such  an  assumption  may  be  subject  to  criticism.  However, 
to  be  consistent  with  current  practice  and  to  use  available  POO  curves, 
the  assumption  will  be  employed  first  in  the  following  formulation.  The 
mathematical  formulation  in  which  the  results  of  multiple  inspections 
are  not  statistically  independent  is  given  in  Appendix  A.  The  issue  of 
dependent  and  independent  multiple  inspections  along  with  the  application 
to  gas  turbine  engine  components  will  be  discussed  later. 

The  formulation  and  solution  will  be  established  in  an  appropriate 
perspective  for  multiple  inspection  procedures.  With  the  mathematical 
solution  put  into  an  appropriate  setting,  one  can  manipulate  various  NDE 
capabilities  (or  POO  curves)  to  achieve  a  most  economical  multiple 
inspection  system  in  terms  of  Type  I  error.  Type  II  error  and  both. 

Let  the  following  POO  curves  associated  with  m  inspection  systems 
be  given; 

P0D(a;l)  =  probability  of  detecting  the  crack  length,  a,  under 
No.  1  NOE  system 

P0D(a;j)  =  probability  of  detecting  the  crack  length,  a,  under 
No.  j  NDE  system;  j=l,2,...,m. 

In  general,  P0D(a;i)  is  different  from  P0D(a;j).  However,  they  may 
also  be  identical  as  a  special  case.  Two  basic  inspection  rules  and 
their  combinations  are  described  in  the  following. 
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1.  UNION  RULE 

A  component  (or  a  critical  location)  is  rejected  if  it  is  rejected 
by  either  one  of  the  NDE  systems  employed  as  shown  in  Figure  4. 

Let  P0D(a;lU2)  =  probability  that  the  crack  length,  a,  is  detected 
by  either  No.  1  NDE  system  or  No.  2  NDE  system  or  both.  Then,  based  on 
the  assumption  that  the  inspections  by  both  NDE  systems  are  statistically 
independent,  the  resulting  POD  curve  for  the  component  is  given  by 

P0D(a;lU2)  =  1  -  [l-P0D(a;l)][l-P0D(a;2)]  (1) 

or 

P0D(a;lU2)  =  P0D(a;l)  +  P0D{a;2)  -  P0D(a;l)P0D(a;2)  (2) 

In  a  similar  fashion,  the  probability  of  detection  (POD)  under  m  NDE 

m 

systems,  denoted  by  P(a;  U  j)  is  obtained  as 

j=l 


m  m 

P0D(a;  U  j)  =  1  -  n  [l-P0D(a;j)]  (3) 

j=l  j*l 


In  general,  the  application  of  multiple  inspections  using  the  union 
rule  will  reduce  the  Type  I  error  and  the  inspection  limit  a^jQ^,  thus 
improving  the  structural  safety  and  reliability.  However,  it  may  also 
increase  the  Type  II  error;  the  extent  of  which  depends  on  the  particular 
functional  form  of  P0D(a;j);  j=l,2,...,m.  An  exceptional  case  is  that 
if  P0D(a;2)  has  a  lower  bound  at  i.e.,  P0D(a;2)  is  equal  to  zero 

for  a^aj^Qg,  then  there  is  no  increase  in  the  Type  II  error  under  two 
inspections.  A  schematical  flow  chart  for  the  inspection  procedures 


using  two  NDE  systems  and  the  union  rule  is  shown  in  Figure  5. 


POD(a;l)  POD(a;2) 

INTERSECTION  RULE 


Figure  4.  Rejected  Components  Shown  in  Shaded  Areas. 
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Figure  5.  Inspection  Procedure  for  Two  NDE  Systems  with  Union  Rule 
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2.  INTERSECTION  RULE 

A  component  (or  a  critical  location)  is  rejected  if  and  only  if  it 
is  rejected  by  all  the  NOE  systems  employed,  as  shown  in  Figure  4. 

If  P0D(a;ln2)  denotes  the  probability  of  detecting  the  crack  length 
a,  by  both  No.  1  and  No.  2  NOE  systems,  then  we  have 

P0D(a;lfi2)  =  P0D{a;l)P0D(a;2)  (4 

In  a  similar  manner,  the  probability  of  detecting  the  crack  length,  a, 
by  every  one  of  m  NOE  systems  employed  is  given  by 


P0D(a;  n^j) 
J=1 


n  P0D(a;j) 
j=l 


In  general,  the  application  of  multiple  inspections  using  the 
intersection  rule  alone  will  degrade  the  NOE  capability.  It  is 
precisely  due  to  such  a  property  that  the  Type  II  error  can  be  reduced. 
However,  caution  should  be  taken  such  that  the  degradation  for  the  Type 
I  error  will  be  insignificant.  For  instance,  using  two  inspection 
systems,  the  condition  for  not  having  a  serious  adverse  effect  on  the 
Type  I  error  is  that  P00(a;2)  should  be  very  close  to  unity  at  as 

will  be  described  later.  If  P0D(a;2)  has  an  upper  bound  at  aj^p^,  i.e., 
P0D(a;2)=l  for  a^a^pg,  then  the  Type  I  error  will  not  be  effected.  A 
flow  chart  for  the  inspection  procedure  using  two  NDE  systems  and  the 
intersection  rule  is  shown  in  Figure  6. 

3.  COMBINATION  RULE  FOR  THREE  INSPECTIONS 

As  described  previously,  the  application  of  either  the  union  rule 
or  the  intersection  rule  alone  cannot  reduce  both  Type  I  and  Type  II 
errors  simultaneously.  However,  a  cOTibined  use  of  both  union  and  inter¬ 
section  rules,  along  with  an  appropriate  choice  of  P0D(a;3),  can  result 
in  a  simultaneous  improvement  for  both  types  of  error. 

Because  of  practical  limitations,  such  as  the  facility  and 
inspection  cost,  we  shall  describe  only  a  possible  combination  of  union 
and  intersection  rules  using  three  inspections  as  follows. 


NO.  1  NDE  SYSTEM  RGECTED  _ A  NO.  2  NDE  SYSTEM 


Figure  6.  Inspection  Procedure  for  Two  NDE  Systems  with 
Intersection  Rule. 
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Let  POD  [a; (lU2)a33  be  the  probability  of  detecting  the  crack  length, 
a,  under  three  inspections,  where  the  union  rule  is  applied  to  No.  1  and 
No.  2  NDE  systems  and  the  intersection  rule  is  employed  for  the  No.  3 
NOE  system.  The  resulting  POD  curve  from  the  application  of  No.  1  and 
No.  2  NDE  systems,  P00(a;lU2),  is  given  by  Equation  2.  Hence,  it  follows 
from  Equation  4  that 

P0D[a;(lU2)n3]  =  P0D(a;3)P0D{a;lU2)  (6) 

Substituting  Equation  2  into  Equation  6,  one  obtains 

P0D[a;(lU2)D3]  =  P0D(a;3) [P0D(a;l)  +  P0D(a;2) 

-  P0D(a;l)P00(a;2)]  (7) 

The  inspection  procedure  for  three  NDE  systems,  presented  in 
Equation  7,  is  shown  by  a  flow  chart  in  Figure  7  and  explained  as  follows 
After  the  first  inspection  by  the  No.  1  NDE  system,  components  are 
divided  into  two  populations:  an  accented  one  and  a  rejected  one.  The 
accepted  population  is  further  inspected  by  the  No.  2  NDE  system  (to 
reduce  Type  I  error)  and  the  accepted  components  are  returned  to  service. 
The  components  rejected  by  the  No.  2  NDE  system  along  with  the  components 
rejected  by  the  No.  1  NDE  system  are  then  inspected  by  the  No.  3  NDE 
system  (to  reduce  Type  II  error).  Then  the  rejected  components  (by  No.  3 
NDE  system)  are  replaced  or  repaired  and  the  accepted  components  are 
returned  to  service.  It  should  be  mentioned  that  P0D(a;3)  for  the  third 
NDE  system  should  be  very  close  to  1.0  (such  as  0.99)  at  9=3^^^,  so  that 
the  adverse  effect  on  the  Type  I  error  is  minimal.  Likewise,  the  band¬ 
width  of  P0D(a;3)  should  be  as  narrow  as  possible  to  reduce  the  Type  II 
error  effectively.  Inspections  using  such  an  NDE  system  may  be  expensive 
and  time  consuming,  but  the  number  of  components  to  be  inspected  by  the 
No.  3  NDE  system  may  be  small.  This  will  be  discussed  later. 

If  m  NDE  systems  are  employed,  then  there  are  many  different 
combinations  of  union  and  intersection  rules  that  can  be  investigated. 
However,  the  basic  mathematical  approach  to  derive  the  resulting  POD 
curve  is  the  same  as  that  described  previously.  For  instance,  for  a 
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special  case  in  which  P0D(a;l)=P0D(a;2)=P0D(a;3),  then  the  resulting  POD 
curve  using  the  union  rule  is  given  by 


P0D(a;  U  j)  =  1  -  [l-PODla;!)]"^  (8) 

j=l 

When  the  fourth  NOE  system  with  P0D(a;4)  is  further  applied  to  inspect 
those  rejected  components  using  the  intersection  rule,  the  resulting  POD 
curve  is  given  by 


3  3 

P0D[a;(  U  j)5l4]  =  P0D(a;4)P0D{a;  U  j) 

j=l 

3 

in  which  P0D(a;  U  j)  is  given  by  Equation  8. 
j=l 

4.  INSPECTION  SEQUENCE  FOR  MINIMUM  NUMBER  AND  COST  OF  INSPECTIONS 


(9) 


The  resulting  POD  curves  under  multiple  inspections  derived  in  the 
previous  sections  are  independent  of  the  sequence  (or  order)  of  appli¬ 
cations  of  multiple  NOE  systems.  However,  the  number  of  inspections 
required  for  components  and  hence  the  cost  of  inspection  is  indeed 
influenced  by  the  sequence  of  inspections.  For  the  case  of  two  in¬ 
spections  with  the  union  rule  as  shown  in  Figure  5,  P0D(a;l)  should  be 
better  than  P0D(a;2)  such  that  fewer  components  will  be  accepted  by  the 
No.  1  NOE  system.  Therefore,  the  number  of  components  to  be  inspected 
by  the  No.  2  NOE  system  is  minimal.  In  general,  for  multiple  inspections 
with  the  union  rule,  the  NOE  system  with  the  highest  resolution  capa¬ 
bility  should  be  applied  first  in  order  to  minimize  the  subsequent  number 
of  inspections.  However,  the  minimization  of  the  number  of  inspections 
implies  the  minimization  of  the  inspection  cost  only  when  the  cost  per 
inspection  for  both  No.  1  and  No.  2  NOE  system  is  almost  the  same.  If 
there  is  a  difference  in  the  inspection  cost  for  two  NOE  systems,  then 
the  system  with  the  lower  inspection  cost  should  be  the  No.  1  NDE  system. 
This  is  because  the  No.  1  NDE  system  has  to  inspect  all  the  components 
whereas  only  those  components  accepted  by  the  No.  1  NDE  system  will  be 
inspected  by  the  No.  2  NDE  system. 


On  the  other  hand,  for  the  case  of  two  inspections  using  the  inter¬ 
section  rule,  as  shown  in  Figure  6,  P0D(a;2)  should  be  better  than  P0D(a;l), 


AFWAL-TR-82-4111 


so  that  the  number  of  components  rejected  by  the  No.  1  NDE  system,  which 
should  be  inspected  by  the  No.  2  NDE  system,  will  be  minimal.  For  the 
case  of  three  inspections  using  the  union-intersection  rule,  the  number 
of  components  to  be  inspected  by  the  No.  3  NDE  system  is  independent  of 
the  inspection  sequence.  However,  for  the  first  two  inspections  using 
the  union  rule,  P0D{a;l)  should  be  better  than  P0D{a;2)  to  achieve  a 
minimum  number  of  inspections. 

Again,  depending  on  the  cost  of  inspection  for  each  NDE  system,  the 
rationale  described  previously  can  be  applied  to  minimize  the  inspection 
cost. 

5.  CORRELATED  MULTIPLE  INSPECTIONS 

The  solutions  obtained  previously  for  the  resulting  POD  curves  under 
multiple  inspection  procedures  are  based  on  the  assumption  that  inspection 
results  from  multiple  NDE  systems  are  statistically  independent.  This 
assumption  is  consistent  with  the  current  practice  for  establishing  the 
POD  curve  for  each  NDE  system  in  which  independent  sampling  has  been 
assumed,  i.e.,  each  inspection  result  from  the  same  NDE  system  is  assumed 
to  be  statistically  independent. 

It  is  more  economical  to  perform  multiple  inspections  using  the 
same  NDE  system.  However,  the  question  of  whether  the  results  of  multiple 
inspections  using  the  same  NDE  system,  but  under  different  inspection 
conditions,  will  be  independent  or  correlated  has  not  yet  been  fully 
resolved.  Take  eddy-current  inspection  for  instance.  The  following 
conditions  may  be  different  in  multiple  inspections;  inspector,  gain  of 
NDE  signal,  scanning  speed  for  inspection,  position  of  probe,  signal 
data  processing,  surface  preparation  of  inspected  parts,  loading  condition 
of  parts,  the  same  system  at  different  locations  or  facilities,  etc. 

For  a  fully  automated  NDE  system  currently  under  development  for 
the  RFC  system,  it  is  anticipated  that  the  error  or  uncertainty  due  to 
human  operation  and  others  will  be  greatly  reduced,  and  the  systematic 
error  or  uncertainty  of  the  NDE  system  itself  will  prevail.  Thus, 
multiple  inspections  using  such  a  fully  automated  NDE  system  alone  for  a 
single  location  may  be  highly  correlated  if  the  preparation  for  the  sur¬ 
face  condition  of  the  location  is  identical. 
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The  mathematical  solutions  for  correlated  inspections  are  presented 
in  Appendix  A.  With  correlated  inspection  systems,  however,  additional 
POD  information  is  needed.  For  instance,  the  conditional  probability  of 
detecting  the  crack  length,  a,  by  the  No.  2  NDE  system  under  the  condition 
that  the  crack  has  been  detected  by  the  No.  1  NDE  system  is  required. 

The  procedure  to  establish  such  a  conditional  POD  curve  is  described  in 
Appendix  B. 

6.  APPLICATION  OF  INDEPENDENT  INSPECTIONS 

One  of  the  motivations  of  the  present  investigations  is  the  appli¬ 
cation  to  the  Retirement-For-Cause  of  gas  turbine  engine  components, 
such  as  disks.  A  disk  usually  contains  many  holes,  such  as  bolt  holes, 
cooling  holes  and  rim  holes,  in  which  cracks  may  occur.  Since  the  crack 
in  each  hole  may  have  a  different  length,  orientation  and  geometry,  the 
inspection  of  one  hole  can  be  assumed  to  be  statistically  independent  of 
the  inspection  of  another  hole.  Likewise,  in  order  to  reduce  Type  II 
error,  the  rejected  holes  may  be  cleaned,  polished  or  even  replicated, 
in  which  case  a  high  resolution  capability  for  the  POD  curve  can  be 
achieved.  It  may  be  reasonable  to  assume  that  the  inspection  for  a  hole 
with  replication  is  statistically  independent  of  the  inspection  for  the 
same  hole  without  replication  even  when  the  same  NDE  system  is  used. 

The  application  of  the  present  results  to  RFC  of  engine  components  will 
be  presented  in  later  sections. 

7.  TYPE  I  AND  TYPE  II  ERRORS 

The  effect  of  multiple  inspections  on  the  POD  curve.  Type  I  and 
Type  II  errors,  and  the  inspection  limit  aj^p^  will  be  demonstrated  later. 
While  the  capability  of  a  particular  NDE  system  is  defined  by  its  POD 
curve,  the  Type  I  and  Type  II  errors  depend  not  only  on  the  POD  curve 
itself  but  also  on  the  pre-inspection  distribution  of  the  flaw  length  in 
the  component.  For  instance,  if  all  the  crack  lengths  in  the  component 
prior  to  inspection  are  smaller  than  then  the  Type  I  error  is  zero. 

Let  F(a)  and  f(a)  be  the  distribution  function  and  the  probability 
density  function,  respectively,  of  the  flaw  length  in  the  component  prior 
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to  inspection.  Then  the  Type  I  and  Type  II  errors,  denoted  by  Pj  and 
Pjj,  respectively,  are  given  by 

CO 

Pj=/  f(x)[l-POD{x)]dx  (10) 

^NDE 

®NDE 

PlI  =/  f(x)P0D(x)dx  (11) 

0 

in  which 

Pj  =  probability  of  missing  (or  accepting)  a  crack  length  longer 
than 

Pjj  =  probability  of  detecting  (or  rejecting)  a  crack  length 
smaller  than 

and  POD(a)  is  the  POD  curve  of  a  particular  NDE  system. 

Both  Pj  and  Pjj,  given  by  Equations  10  and  11,  are  the  quantitative 
measures  of  Type  I  and  Type  II  errors.  Two  qualitative  measures  of  Type 
I  and  Type  II  errors  which  depend  exclusively  on  the  POD  curve,  may  also 
be  appealing. 


Aj  =/  [l-P0D(x)]dx  (12) 

®NDE 
^NDE 

Ajj  */  P0D(x)dx  (13) 

0 

It  is  apparent  from  Equations  12  and  13  that  Aj  is  the  area  above  the  POD 
curve  from  aj^p^  to  infinity,  and  Ajj  is  the  area  under  the  POD  curve  from 
zero  to  as  shown  in  the  shaded  areas  of  Figure  3.  While  Aj  and  Ajj 

are  not  the  quantitative  measures  of  Type  I  and  Type  II  errors,  they  may 
serve  for  qualitative  comparisons  between  the  capability  of  various  NDE 
systems  when  the  value  of  the  inspection  limit,  a^^p^,  is  fixed. 


When  multiple  inspections  are  employed,  the  PQD(a)  function 
appearing  in  Equations  10-13  represents  the  resulting  POD  curve  derived 


in  Equations  1-9.  Hence,  it  should  be  replaced  by  the  approi 
corresponding  POD  curve  resulting  from  multiple  inspections. 
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8.  DISTRIBUTION  OF  PRE-INSPECTION  FLAW  LENGTH 


Theoretically,  the  distribution  of  the  flaw  length  for  a  component 
prior  to  inspection  can  be  derived  from  the  distributions  of  times  to 
crack  initiation  (or  equivalent  initial  flaw  length),  crack  growth  rates 
in  service,  and  service  loads.  Such  an  approach,  however,  is  very 
complex.  From  the  NOE  standpoint,  it  is  reasonable  to  assume  the  Weibull 
distribution  for  pre-inspection  flaw  lengths  for  illustration.  Information 
on  the  statistical  variabilities  of  the  equivalent  initial  flaw  length, 
the  crack  propagation  rate  and  the  service  loads  for  engine  disks 
indicates  that  the  statistical  dispersion  (coefficient  of  variation)  of 
the  crack  length  in  service  is  of  the  order  of  100%.  The  Weibull  distri¬ 
bution  with  a  100%  coefficient  of  variation  degenerates  into  a  special 
case  of  the  negative  exponential  distribution.  As  a  result,  the  distri¬ 
bution  of  the  pre-inspection  flaw  length  may  be  assumed  to  follow  the 
negative  exponential  distribution,  i.e.. 


f(a)  =  X  ;  a^O 

F(a)  =  1  -  e"^^  ;  a^O 


(14) 


in  which  f(a)  and  F(a)  are,  respectively,  the  probability  density  function 
and  the  distribution  function  of  the  pre-inspection  flaw  length.  In 
Equation  14,  l/x  represents  the  average  flaw  length. 


The  lognormal  distribution  may  also  be  used  for  the  pre-inspection 
flaw  length  distribution  as  follows; 

2 

a  -  M 


f(a)  .  -i°a-g_expi  -  l(  ’°9  ^  -  ") 

/ST  aa  ^  \  ”  /  \ 


(15) 


F(a)  =  ^ 


in  which  p  and  a  are  the  mean  value  and  the  standard  deviation, 
respectively,  of  the  logarithm  of  the  flaw  length,  and  $(  )  is  the 
standardized  normal  distribution  function. 
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SECTION  III 
NUMERICAL  EXAMPLES 


1.  EXAMPLE  NO.  1 

Eddy  current  inspection  data  for  fastener  holes  in  skin  and  stringer 
wing  assembly,  referred  to  as  HAVE  CRACK,  were  available  in  Reference  3 
and  analyzed  by  the  binomial  method.  The  same  data  set  was  further 
analyzed  in  Reference  8  using  the  regression  method  with  the  assumed 
functional  form  for  the  POD  curve  as  follows; 

P0D(a;l}  =  exp(a  Btna) 

1  +  exp(a  +  atna) 

in  which  a=-2.9  and  e=1.7  (Figure  8).  This  POD  curve  is  replotted  in 
Figure  9(a)  as  Curve  1.  The  crack  length  associated  with  90X  and  96X 
detection  probabilities  are,  respectively,  20.05  mm  and  35.7  ram. 

The  data  set  used  to  establish  the  POO  curve  was  generated  using 
the  same  cracked  specimens  but  inspected  by  different  inspectors  and  NDE 
systems  at  different  locations.  If  the  components  are  inspected  twice 
at  different  locations,  i.e.,  P00(a;l)=P0D(a;2),  the  resulting  POD  curve, 
i.e.,  P0D(a;lU2),  with  the  union  rule  is  shown  as  Curve  2  (Equation  2) 
in  Figure  9(a).  One  can  further  improve  the  Type  I  error  by  performing 
a  third  inspection  using  the  union  rule.  With  P0D(a;l)=P0D(a;2)®P0D(a;3), 
the  resulting  P0D(a;lU2U3) ,  Equation  3,  is  displayed  in  Figure  9(a)  as 
Curve  3.  If  the  inspection  limit,  is  required  to  be  12.4  tnm,  i.e,, 

aNDE=12.4  mm,  then,  Aj  and  Ajj  (Equations  12  and  13)  are  computed  and 
shown  in  Table  1.  It  is  observed  from  Table  1  that  Aj  reduces  and  Ajj 
increases  as  the  number  of  inspections  increases. 

The  distribution  of  the  flaw  length  prior  to  inspection  is  given  by 
Equation  14  in  which  the  average  flaw  length  is  assumed  to  be  5  mm,  i.e., 

X  =  0.2/mm.  The  probability  density  function  given  by  Equation  14  is 
displayed  in  Figure  9(a)  as  a  dashed  curve.  With  such  a  pre-inspection 
flaw  length  distribution,  the  probability  that  the  crack  length  will 
exceed  a|n[)£®12.4  mm  is  8.37%,  i.e.,  on  the  average  there  are  8,37%  of  the 
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Figure  8.  Eddy  Current  Inspections  of  Skin  and  Stringer  Wing 
Assembly,  60  Inspections  Per  Fastener  Hole 
(from  References  3  and  8). 
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components  that  will  have  a  crack  length  longer  than  12.4  mm. 

Therefore,  the  average  percentage  of  good  components,  Pg,  prior  to 
inspection  is  91.63X,  i.e.,  Pg=91.36X. 

The  Type  I  and  Type  II  errors,  Pj  and  Pjj,  are  computed  from 
Equations  10  and  11  in  which  POD(x)  is  replaced  by  P0D(x;l),  P0D(x;lU2) 
and  POD{x;lU2U3),  respectively.  The  results  are  shown  in  Table  1.  It 
is  observed  from  Table  1  that  the  Type  I  error  reduces  drastically  as 
the  number  of  inspections  with  the  union  rule  increases.  However,  the 
Type  II  error  increases  simultaneously. 


TABLE  1 

TYPE  I  AND  TYPE  II  ERRORS;  a^pg  =  12.4  mm,  Pg  =  91.63X 

TABLE  1 


No.  of  Inspect! 

ons 

*1 

^I 

Pl 

PlI 

Type 

(mm) 

(mm) 

1 

P0D(a;l) 

2.5 

5.99 

1.15x10"^ 

28.6% 

2 

P0D(a;lU2) 

0.24 

8.30 

1 .72x10’^ 

42.7% 

3 

P0D[a;(lU2)n3*] 

0.24 

1.1 

1.72x10"^ 

5.9% 

3 

P0D(a;lU2U3) 

0.03 

9.36 

2.74x10'^ 

50.8% 

4 

P0D[a;(lU2U3)n3 

*] 

0.03 

1.38 

2.74x10"^ 

6.2% 

HE 


E 
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additional  inspection  using  the  intersection  rule  is  assumed  to  resemble 
that  of  the  x-ray  system  shown  in  Figure  1,  i.e.. 


P0D(a;3*)  =  0 

_  a-9.6 

2.8 


for  a<9.6  mm 


for  9.6<a<l2.4 


for  a>12.4 


It  is  obvious  that  the  POD  curve  given  above  is  a  straight  line  between 
9.6  mm  and  12.4  mm  (Figure  1). 

The  resulting  POD  curve,  denoted  by  P0D[a;(lU2)n3*]  and  P0D[a; (11)203) 
Ji3*],  presented  in  Figure  9(b)  as  solid  and  dashed  curves,  respectively. 
Note  that  P00[a; (102)^3*]  is  the  resulting  POD  curve  under  two  inspections 
with  the  union  rule,  i.e.,  P0D(a;l )=P0D(a;2) ,  as  well  as  the  third  inspection 
with  the  POD  curve  given  by  Equation  17.  Lilcewise,  P0D[a; (lU2U3)fi3*]  is 
the  resulting  POD  curve  under  three  inspections  with  the  union  rule,  and 
an  additional  inspection  with  the  intersection  rule,  where  the  POD  curve 
is  given  by  Equation  17.  Both  P0D[a;(lU2)O3*]  and  P0D[a;(lU2U3)o3*]  are 
computed  from  Equations  7  and  9,  respectively. 

The  Type  I  and  Type  II  errors  (Pj  and  Pjj),  as  well  as  Aj  and  Ajj, 
are  computed  from  Equations  10-13  and  the  results  are  shown  in  Table  1. 

It  is  observed  from  Table  1  and  Figure  9(b)  that  multiple  inspections 
can  reduce  both  Type  I  and  Type  II  errors  simultaneously. 

2.  EXAMPLE  NO.  2 

Eddy  current  inspections  of  etched  fatigue  cracks  in  2219-T87 
aluminum  plates  are  available  from  Reference  2,  in  which  the  optimized 
probability  method  along  with  the  binomial  approach  was  used  to  establish 
the  95%  confidence  POD  curve  shown  in  Figure  10.  For  situations  in  which 
both  Type  I  and  Type  II  errors  are  important,  such  as  the  Retirement- 
For-Cause  of  engine  components,  the  mean  POD  curve  should  be  used.  For 
the  purpose  of  illustration,  the  mean  POO  curve  is  approximated  by 


P0D(a;l)  =  1-e"'^  ;  a>0 


PROBABILITY  OF  DETECTION(%) 


CRACK  LENGTH(mm) 

Figure  10.  Optimized  Probability  Method  -  361  Eddy  Current 
Inspections  of  Etched  Fatigue  Cracks  in  2219-T87 
Aluminum  Plates  (from  Reference  2). 

with  3^=0. 46/mm  as  shown  in  Figure  11  (Curve  1).  From  Curve  1  it  is 
observed  that  the  inspection  limit  corresponding  to  a  90%  detection 
probability  is  5  mm,  i.e.,  5  mm.  In  what  follows,  5  mm  will  be 

used  as  the  inspection  limit  or  the  reset  crack  length, 

Suppose  the  POD  curve  for  the  second  NDE  system  is  given  by 


P0D(a;2)  =  1  -  e 
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with  62*0* P0D(a;2)  given  by  Equation  19  is  plotted  in 
Figure  11  as  Curve  2.  It  is  obvious  that  the  second  NOE  system  has  a 
worse  detection  capability  than  the  first  one. 


CRACK  LENGTH,  MM 

Figure  11.  POD  Curves  (1  and  2)  and  Resulting  POD  Curves  Using 
the  Union  Rule. 

With  the  application  of  the  second  NDE  system  using  the  union  rule, 
the  resulting  POD  curve,  i.e.,  P0D(a;lU2),  is  shown  in  Figure  11  as  Curve 
3.  Now  the  inspection  limit,  aj^oE"^  resulting  POD  curve 

corresponds  to  a  96<  detection  probability.  Thus,  it  is  shown  that  even 
with  a  poor  NDE  system,  a  second  inspection  using  the  union  rule  can 
reduce  the  Type  I  error.  The  measures  Aj  and  Ajj  for  Type  I  and  Type  II 
errors  defined  in  Equations  12  and  13  are  obtained  as  follows; 

Aj  =  6xp  [-(bj  ^2^^NDE^^^^1  ^2^  (20) 

^II  ~  ^NDE  ”  ^  S2^*NDE^ ^^^®1  ^  ^2^ 

Assuming  that  the  second  NDE  system  has  the  same  resolution  capa¬ 
bility  as  the  first  one,  i.e.,  s resulting  POD  curve, 
i.e.,  P0D(a;lU2),  is  plotted  in  Figure  11  as  Curve  4.  Now  the  detection 
probability  at  the  inspection  limit  aMr,r=5  mm  increases  to  99%. 


DENSITY  FUNCTION 
OF  PREINSPECTION  FLAWS 
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The  distribution  of  the  pre-inspection  flaw  length  is  assumed  to 
follow  the  negative  exponential  distribution  given  by  Equation  14 
with  x=0.7/mm  which  is  also  shown  in  Figure  11  (dashed  curve).  The 
probability  that  a  crack  length  will  exceed  a|^o£=5  mm  is  computed  as  3515. 
In  other  words,  on  the  average,  3515  of  the  components  should  be  repaired 
or  replaced  if  an  ideal  NDE  system  is  used,  i.e.,  Pq=97%.  The  Type  I 
and  Type  II  errors  given  in  Equations  10  and  11  are  obtained  as 


^I  "  A+e^+e2  j 

^II  "  ]  ■  X+6^+32  ‘  jj  (23) 


Values  of  Aj,  Ajj,  Pj  and  Pjj  are  shown  in  Tables  2  and  3. 

As  expected,  the  Type  I  error  is  reduced  under  two  inspections  using 
the  union  rule  whereas  the  Type  II  error  increases.  From  Tables  2  and 
3,  the  Type  II  error,  either  under  single  or  double  inspections,  may  not 
be  acceptable  from  an  economical  standpoint.  As  a  result,  a  third  inspec 
tion  is  applied  using  the  intersection  rule  as  follows. 


The  third  NDE  system  should  have  a  high  detection  probability  at 
the  inspection  limit,  af^Q£=5  mm,  and  it  should  also  be  narrow-banded  to 
be  effective  for  reducing  the  Type  II  error.  Such  a  high  quality  POD 
curve  may  be  achieved  using  the  first  NDE  system  if  the  cracked  locations 
are  cleaned,  polished  or  replicated.  For  illustrative  purposes,  the  POD 
curve,  P0D(a;3),  is  assumed  to  be  of  the  following  forms; 


P0D(a;3)  =  0 


a<ai 


a-a. 


a^-ai 


a^<a^a2 


a>ar 


or 


P0D(a;3)  =  1  -  exp  ]  -(a/6f.) 


"0 


a^O 


(24) 


(25) 
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TABLE  2 

TYPE  I  AND  TYPE  II  ERRORS;  =  0.46/nm,  3^  “  0.198/mm, 
^NOE  =  5  mm,  Pg  =  97X 


No.  of 

Inspections 

Rule 

(mm) 

(mm) 

Pi 

Pii 

1 

3.04 

1 .83x10"^ 

36% 

2 

U 

0.066 

3.54 

0.57x10"^ 

45% 

★ 

3 

U-1 

0.066 

1.16 

4.7% 

3** 

U-1 

0.069 

1 .04 

3.84% 

*:  POO(a;3)  given  by  Equation  24  with  ai=2.5  mm,  a2=5.0  mm 


**;  P0D(a;3)  given  by  Equation  25  with  aQ=8.0,  3q=4.13  mm 


TABLE  3 


TYPE  I  AND  TYPE  II  ERRORS;  =  0.46/mm,  B^  =  0.46/mm, 
®NDE  =  5  mm,  pQ  =  97% 


No.  of 

Inspections 

Rule 

*1 

(mm) 

^I 

(mm) 

Pi 

PlI 

1 

0.218 

1 .83x10"^ 

36% 

2 

U 

0.011 

3.92 

0.13x10"^ 

54% 

★ 

3 

U-1 

0.011 

1 .22 

0.13x10"^ 

5% 

3 

U-I 

0.0112 

1.09 

0.15x10"^ 

4% 

*:  P0D(a;3)  given  by  Equation  24  with  a^=2.5  mm,  a2=5.0  mm 
**:  P0D(a;3)  given  by  Equation  25  with  0^=8. 0,  3q=4.13  mm 
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The  POO  curve,  P0D(a;3),  given  by  Equation  25  is  a  Weibull  function 
whereas  the  curve  given  by  P0D(a;3)  in  Equation  24  is  a  straight  line 
between  a^  and  ^2  which  resembles  the  curve  for  the  x-ray  system  shown 
in  Figure  1.  Both  Equations  25  and  24  are  plotted  in  Figure  12  as  Curve 
1  and  Curve  2,  respectively,  with  oiq=8  and  Bq=4.13  mm,  a2^=2.5  mm  and 
a«=5.0  mm. 


CRACK  LENGTH,  MM 

Figure  12.  POD  Curves  for  the  Third  Inspection  System  P0D{a;3). 


The  resulting  POD  curves,  i.e.,  P0D[a; (lU2)n3]  given  by  Equation  9, 
are  plotted  in  Figures  13(a)-13(d)  for  different  combinations  of  P0D(a;1), 
P0D(a;2)  and  P0D(a;3),  given  by  Equations  18,  19,  24  and  25.  These  figures 
demonstrate  a  significant  improvemert  of  the  POD  curve  using  three  inspec¬ 
tions  with  the  union-intersection  rule  described  previously.  The  corre¬ 
sponding  Aj,  Ajj,  Pj  and  Pjj  values  are  shown  in  Tables  2  and  3.  As 
observed  from  these  tables,  both  Type  I  and  Type  II  errors  have  been 
reduced  significantly  through  multiple  inspections.  It  is  further  noticed 
that  the  values  of  Aj  and  ^II  serve  only  as  qualitative  measures  of  Type 
I  and  Type  II  errors,  since  they  are  not  directly  proportional  to  Pj  and 
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CRACK  LENGTH,  MM 

Figure  13(a).  Resulting  POD  Curves  Under  Three  Inspections  Using  the 
Union-Intersection  Rule:  (a)  P0D[a; (lU2)o3], 

Weibull  P0D(a;3). 


Figure  13(b). 


CRACK  LENGTH,  MM 

Resulting  POD  Curves  Under  Three  Inspections  Using  the 
Union-Intersection  Rule:  (b)  P0D[a;(lU2)n3] , 

Straight  Line  P00(a;3). 
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5  6 

CRACK  LENGTH,  MM 


Figure  13(c).  Resulting  POD  Curves  Under  Three  Inspections  Using  the 
Union-Intersection  Rule:  (c)  P0D[a; (lUl)ft3], 

Weibull  P0D(a;3). 


0  I 


Figure  13(d) 


CRACK  LENGTH,  MM 

Resulting  POD  Curves  Under  Three  Inspections  Using  the 
Union-Intersection  Rule:  (d)  P0D[a; (lUl)n3], 

Straight  Line  P0D(a;3). 
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Tables  2  and  3  indicate  a  slight  increase  in  the  Type  I  error  when 
the  Weibull  POD  function  is  applied,  (Equation  25)  for  the  third 
inspection.  This  is  because  the  detection  probability  at  the  inspection 
limit,  a|^Q£,  in  Equation  25  is  slightly  less  than  1.0. 

Because  of  the  reduction  in  Type  I  error  using  the  union  rule  for 
P0D(a;l)  and  P0D(a;2)  it  is  possible  to  reduce  the  inspection  limit 
^NDE*  Suppose  the  inspection  limit  is  reduced  to  3.5  mm  which  corresponds 
to  a  96%  detection  probability  as  shown  by  Curve  4  of  Figure  11.  Then 
the  average  percentage  of  good  components  is  Pq=91.4%.  The  P0D(a;3) 
function  should  be  chosen  such  that  it  is  close  to  unity  at  a=3.5  mm. 
Again,  Equations  25  and  24  are  used  with  0^=8. 0  and  6q=2.89  nm,  aj=1.0 
mm,  and  32=3.5  mm.  These  POD  curves  are  shown  in  Figure  12  as  Curve  3 
and  Curve  4,  respectively.  The  results  of  Aj,  Ajj,  Pj  and  Pjj  are 
presented  in  Table  4.  As  indicated  by  Table  4,  multiple  inspection 
procedures  are  capable  of  reducing  both  Type  I  and  Type  II  errors. 


TABLE  4 

TYPE  I  AND  TYPE  II  ERRORS;  Sj  =  0.46/mm,  $3  =  0.46/mm, 
a^rjr  =  3.5  mm,  P«  =  91.4% 


No.  of 

Inspections 

Rule 

(mm) 

^I 

(mm) 

Pi 

Pii 

1 

0.435 

1.76 

1.0  xio"^ 

32.1% 

2 

U 

0.043 

2.56 

1 .49x10"^ 

48.3% 

* 

3 

U-I 

0.043 

1.124 

12.9% 

** 

3 

m 

0.044 

0.724 

6.37% 

*;  P0D(a;3 

**:  P0D(a;3 

)  given 
)  gi ven 

by  Equation  24  with  3-1=1. 0  mm,  33=3.6  mm 
by  Equation  25  with  aQ=8.0,  3^=2. 89  mm 

3 
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As  a  final  example,  the  POO(a;l)  given  by  Equation  18  is  used  whereas 
the  P0D(a;2)  has  a  lower  bound  at  aQ=1.5  mm,  i.e., 

P0D{a;2)  =  0  for  a<ag 

ft  /  ^  ^  (26) 

1-e’^^  ^or  a>,aQ 

where  0^=0. 46/mm.  Both  P00(a;l)  and  P0D(a;2)  given  by  Equations  18  and 
26  are  plotted  in  Figure  14  as  Curves  1  and  2,  respectively.  The  resulting 
POD  curve  using  the  union  rule,  P0D{a;lU2),  is  also  shown  in  Figure  14 
as  Curve  3  (dashed  curve).  Note  that  P0D(a;lU2)  is  identical  to  P0D(a;l) 
in  the  region  a<aQ=1.5  mm. 

Then,  the  same  P0D(a;3)  given  by  Equations  24  and  25,  respectively, 
are  used  in  conjunction  with  the  intersection  rule.  The  Type  I  and  Type 
II  errors  resulting  from  multiple  inspections  are  shown  in  Table  5. 


Table  5  clearly  indicates  the  advantage  of  multiple  inspections  in  terms 
of  Type  I  and  Type  II  errors. 
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TABLE  5 

TYPE  I  AND  TYPE  II  ERRORS;  =  0.46/mm,  62  =  0.46/mn, 
Sq  —  1*5  irwrij  ^ 


No.  of 

Inspections 

Rule 

^I 

Pii 

1 

0.218 

3.04 

1  .83x10"^ 

36% 

2 

U 

0.022 

3.39 

39.7% 

* 

3 

U-I 

0.022 

1.12 

BlffHTM 

4.6% 

*★ 

3 

U-I 

0.023 

1 .06 

3.9% 

*:  POD(a;3)  given  by  Equation  24  with  a.j=2.5  mm,  a2=5.0  mm 
**:  P0D(a;3)  given  by  Equation  25  with  ctQ=8.0,  eQ=4.13  mm 
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SECTION  IV 

INSPECTION  RELIABILITY  OF  ENGINE  DISKS 

For  gas  turbine  engine  components,  such  as  disks,  inspections  are 
performed  for  each  critical  location,  including  bolt  holes,  cooling  holes, 
rim  holes,  etc.  A  disk  which  consists  of  a  large  number  of  holes  is 
normally  retired  if  at  least  one  hole  is  rejected  by  the  NDE  system. 

The  POD  curve  established  for  any  one  particular  NDE  system  applies  to 
one  hole  only.  Thus,  the  POD  curve  for  a  disk  containing  a  large  number 
of  holes  differs  significantly  from  that  of  the  NDE  system.  As  a  result, 
one  important  task  in  the  Retirement-For-Cause  analysis  is  to  establish, 
if  possible,  the  POD  curve  for  a  disk  from  the  POD  curve  for  a  given  NDE 
system.  An  exploratory  analytical  study  of  this  subject  follows. 

For  simplicity  of  presentation,  the  analysis  is  limited  to  one  type 
of  hole  (such  as  bolt  holes)  to  examine  the  relation  between  the  POD 
curve  of  the  NDE  system  and  that  of  a  disk.  The  general  approach  in 
dealing  with  a  disk  containing  bolt  holes,  cooling  holes,  rim  holes,  and 
other  critical  locations  which  are  subject  to  inspection  will  be  presented 
in  a  future  report. 

The  acceptance  or  rejection  of  a  disk  depends  not  only  on  the  POD 
curve  of  the  NDE  system  employed  but  also  on  the  distribution  of  the 
flaw  length  in  each  hole  prior  to  inspection.  In  this  connection,  two 
extreme  cases  will  be  investigated.  Both  cases  are  mathematically  simple. 

Case  1  -  Completely  Correlated  Crack  Length:  At  one  extreme,  the 
crack  length  in  a  bolt  hole  is  completely  correlated  with  that  in  other 
bolt  holes  in  the  same  disk.  Thus,  the  crack  length  in  each  bolt  hole 
of  a  given  disk  is  identical.  However,  the  crack  length  varies  from 
disk  to  disk  and  it  is  a  random  variable  across  the  ensemble  of  disks  as 
schematically  shown  in  Figure  15(a). 

Case  2  -  Completely  Uncorrelated  Crack  Length;  At  another  extreme, 
the  crack  length  varies  from  one  bolt  hole  to  another  in  the  same  disk. 

The  crack  length  in  each  bolt  hole  is  a  statistically  independent  and 
identically  distributed  random  variable  as  schematically  shown  in  Figure 
15(b). 
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Basically,  the  statistical  variability  of  the  crack  length  in  each 
bolt  hole  is  contributed  to  by  two  major  elements;  (i)  the  statistical 
variability  of  service  loads  and  (ii)  the  statistical  variability  of 
material  properties,  including  the  time  to  crack  initiation,  the  initial 
flaw  length  and  the  crack  propagation.  While  the  service  load  is  a 
statistical  variable,  it  varies  from  engine  (or  disk)  to  engine  (disk) 
and  the  loading  applied  to  each  bolt  hole  in  the  same  disk  is  almost 
identical.  The  material  properties,  however,  vary  from  one  bolt  hole  to 
another  not  only  in  the  same  disk  but  also  amorig  different  disks.  Thus, 
at  one  extreme,  if  the  statistical  variability  of  material  properties  is 
insignificant  compared  to  that  of  the  service  load,  then  the  crack  length 
in  each  bolt  hole  of  the  same  disk  is  almost  completely  correlated  as 
described  in  Case  1.  At  the  other  extreme,  if  the  scatter  of  the  service 
load  is  very  small  compared  to  the  statistical  variability  of  material 
properties,  then  the  crack  length  in  each  hole  will  be  statistically 
independent,  as  described  in  Case  2. 

In  reality,  however,  depending  on  the  mission  characteristics  of  a 
particular  type  of  aircraft,  the  true  situation  will  lie  between  the  two 
extreme  cases  described  previously.  Theoretically,  the  exact  distribution 
of  the  crack  length  in  each  bolt  hole  of  a  disk  prior  to  inspection  can 
be  derived  from  the  statistical  distributions  of  the  time  to  crack 
initiation,  the  equivalent  defect  size,  the  crack  growth  rate  and  the 
service  loads.  However,  such  an  approach  is  analytically  complex.  A 
separate  report  on  an  analytical  RFC  approach  will  be  presented  later. 

From  the  standpoint  of  NOE  reliability  for  engine  components  as  well  as 
the  possible  improvement  of  NDE  capability  through  multiple  inspections, 
these  two  extreme  cases  outlined  herein  are  investigated  because  of  their 
mathematical  simplicity.  Likewise,  a  greater  insight  into  the  influence 
of  NDE  capability  on  the  RFC  of  engine  disks  can  be  achieved  from  the 
present  investigation.  Analytical  solutions  for  the  two  extreme  cases 
are  obtained  in  the  following. 

1.  COMPLETELY  CORRELATED  CRACK  LENGTH 

When  the  crack  length  in  each  bolt  hole  is  identical  in  the  same 
disk,  but  varies  from  disk  to  disk,  the  inspection  of  a  disk  containing 
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m  holes  is  equivalent  to  inspecting  a  crack  length  m  times  (one  time  for 
each  hole)  independently,  using  the  union  rule.  This  is  because  the 
disk  will  be  retired  if  any  hole  is  rejected  by  the  NDE  system. 

Let  P0D(a;l)  be  the  POO  curve  associated  with  the  particular  NDE 
system  employed  for  inspection.  Then  the  POD  curve  for  a  disk  containing 
m  bolt  holes,  denoted  by  P0D(a;Um),  is  obtained  from  P0D(a;l)  using  the 
union  rule. 


P00(a;Um)  =  1  -  [  l-P0D(a;  1 )]  "*  (27) 

A  good  disk  is  defined  as  a  disk  in  which  the  crack  length  in  each 
bolt  hole  is  smaller  than  the  inspection  limit  If  F(a)  and  f{a) 

denote  the  distribution  function  and  the  probability  density  function, 
respectively,  of  the  pre-inspection  flaw  length  in  the  bolt  hole,  then 
the  probability  that  a  disk  is  good  is  given  by 

^NDE 

Pg  =  P[a<a,^Q^]  =J  f(x)dx  =  F(a^p^)  (28) 

0 

in  which  F(3|^0£)  is  the  distribution  function  of  the  pre-inspection  flaw 
length  in  the  hole  evaluated  at  aj^Q^.  Note  that  F(a|^Q^)  can  also  be 
interpreted  as  the  average  percentage  of  good  disks  prior  to  inspection. 

The  Type  I  and  Type  II  errors,  Pj  and  Pjj,  of  a  disk  are  given  by 
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Numerical  examples  will  be  given  to  demonstrate  the  inspection 
reliability  in  terms  of  Type  I  and  Type  II  errors  for  engine  disks. 

Since  Type  I  and  Type  II  errors  depend  on  both  the  POD  curve  of  the  NDE 
system  and  the  distribution  of  the  pre- inspect  ion  flaw  length,  various 
types  of  POD  functions  will  be  investigated.  Thus,  one  can  determine 
the  direction  in  which  the  POD  curve  of  the  NDE  system  should  be  improved 
to  achieve  a  maximum  benefit  from  the  inspection. 

It  should  be  emphasized  that  the  numerical  examples  given  in  Sections 
IV  and  V  are  mainly  for  illustration.  The  range  of  the  flaw  length  for 
POD  curves  employed  is  far  from  the  current  gas  turbine  engine  inspection 
technology,  although  the  POD  curves  are  obtained  from  available  data. 
Furthermore,  the  flaw  length  in  disk  holes  is  not  realistic.  However, 
the  conclusions  and  trends  obtained  for  both  Type  I  and  Type  II  errors 
will  not  be  altered.  This  is  because  when  the  flaw  length  in  holes  and 
the  flaw  length  for  POD  curves  are  reduced  by  the  same  factor  y.  not 
only  the  Type  I  and  Type  II  errors  remain  the  same,  but  also  the  average 
percentage  of  good  holes  at  the  time  of  inspection  maintenance  remains 
identical.  Realistic  examples  for  the  current  state  of  technology  in 
engine  disk  inspections  will  be  presented  in  Section  VI. 

2.  EXAMPLE  NO.  3 

a.  Exponential  POD  Function 

The  exponential -type  POD  curve  associated  with  the  eddy-current 
inspection  of  aluminum  plate  given  by  Equation  18  with  3^=0. 46/mm  is 
considered,  i.e.. 


P0D(a;l)  =  1  -  e  ;  a^O  (30) 

From  Equation  30,  the  detection  probability  for  the  inspection  limit 
^NDE"^  mm  is  90%. 

The  same  pre-inspection  density  function,  f(a),  for  the  flaw 
length  in  each  hole  as  given  by  Equation  14  with  A=0.7/mm  is  used  herein. 
With  such  a  pre-inspection  flaw  length  distribution,  there  are,  on  the 
average,  3%  of  disks  that  will  have  a  crack  length  exceeding  the  inspection 
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limit  amQ£=5  mm,  i.e.,  Pg=97%  in  Equation  28.  The  resulting  POD  curve 
for  the  disk,  P0D(a;U10),  containing  ten  bolt  holes,  i.e.,  m=10  given  by 
Equation  27,  is  plotted  in  Figure  16.  Also  shown  in  Figure  16  are  the 
P0D(a;l)  of  the  NOE  system  and  the  pre-inspection  density  function  of 
the  flaw  length  f(a).  It  is  observed  that  the  bandwidth  of  the  POD  curve 
for  a  disk,  P0D(a;U10),  is  much  smaller  than  that  of  the  POD  curve  of 
the  NDE  system,  P0D(a;l).  With  the  inspection  limit  a|^Q£=5  mm,  the  Type 
I  and  Type  II  errors  (Pj  and  Pjj)»  as  well  as  the  average  percentage  of 
good  disks,  Pg,  are  presented  in  Table  6. 

In  the  first  row  of  Table  6,  the  disk  is  assumed  to  consist  of 

only  one  hole  whereas  the  second  row  represents  the  results  for  ten  bolt 

holes  in  one  disk  (e.g.,  TF  33  disk).  A  comparison  between  the  results 

in  the  first  and  the  second  rows  of  Table  6  indicates  that  the  Type  I 

-3  -13 

error  reduces  drastically  from  1.83x10  for  one  bolt  hole  to  4.1x10 
for  a  disk  with  ten  bolt  holes.  However,  the  Type  II  error  increases 
from  36.8%  to  83%. 

Since  the  Type  I  error  is  extremely  small  for  the  disk,  one  can 
reduce  the  inspection  limit  a^jp^  in  order  to  reduce  the  Type  II  error 
while  keeping  the  Type  I  error  at  an  acceptable  level.  The  results  using 
different  values  of  the  inspection  limit,  are  presented  in  Table 

6.  As  expected,  the  Type  II  error  reduces  and  the  Type  I  error  increases 
as  the  inspection  limit  a|^pg  decreases. 

Note  that  the  Type  II  error  is  the  probability  of  rejecting  a 
disk  that  is  a  good  disk,  i.e., 

Pjj  =  P[RnG]  =  P[R,G]  (31) 

in  which  {R}=the  event  of  rejecting  a  disk  and  {G}=the  event  that  the 

disk  is  good.  As  the  values  of  the  inspection  limit,  a^^p^,  vary,  the 

average  percentage  of  good  disks  changes  accordingly.  Therefore,  another 

measure  of  the  Type  II  error,  referred  to  as  the  conditional  Type  II 

error,  may  be  of  interest.  The  conditional  Type  II  error,  denoted  by 
★ 

Pjj,  is  the  probability  of  rejecting  a  disk  under  the  condition  (or 
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given)  that  the  disk  is  good.  Hence,  we  have  from  Equations  28  and  31 
that 

it 

The  conditional  Type  II  error,  Pj^,  is  also  given  in  Table  6. 

The  results  shown  in  Table  6  are  striking.  It  appears  that  the 
Type  I  error  for  a  disk  is  insignificant  whereas  the  Type  II  error  may 
be  a  serious  problem.  Likewise,  the  inspection  limit  aj^^^  should  be 
determined  by  the  POD  curve  for  the  disk  rather  than  the  POD  curve  of 
the  NDE  system. 

b.  Weibull  POD  Function 

The  sensitivity  of  Type  I  and  Type  II  errors  for  disks  with 
respect  to  different  types  of  POD  function  is  studied  herein.  Let  us 
consider  a  Weibull-type  POD  function  for  the  NDE  system  as  follows; 

“O 

P0D(a;l)  =  1  -  jexp  j  ;  a>0  (33) 

in  which  is  the  shape  parameter  representing  the  bandwidth  of  the  POD 
curve  and  6q  is  the  scale  parameter  representing  the  central  location  of 
the  POD  curve.  To  be  consistent  with  the  previous  example,  we  choose 
oiq=3,0  and  p.q=3.79  mm  such  that  the  inspection  limit,  a|^Q£=5  mm, 
corresponds  to  a  90%  detection  probability.  The  distribution  of  the 
pre-inspection  flaw  length  is  identical  to  that  considered  in  the 
previous  example. 

The  Weibull  POD  curve  (Equation  33),  the  POD  curve  for  the  disk 
P0D(a;U10)  (Equation  27),  and  the  density  function  of  the  pre-inspection 
flaw  length  are  shown  in  Figure  17.  Unlike  the  exponential  POD  function 
given  by  Equation  30,  the  Weibull  function  is  an  S  shaped  function; 
hence  the  areas  Aj  and  Ajj  are  smaller. 

The  results  of  Type  I  and  Type  II  errors,  conditional  Type  II 
error,  and  the  average  percentage  of  good  disks  before  inspection,  are 
shown  in  Table  7  for  different  values  of  the  inspection  limit  The 
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Figure  17.  Weibull  POD  Curves  for  the  NDE  System  and  Disk 

TABLE  7 

TYPE  I  AND  TYPE  II  ERRORS  FOR  DISK:  WEIBULL  POD  CURVE  WITH 
a„  =  3.0.  =  3.79  mm 


No.  of 
Holes 


97.0% 


97.0% 


94.0% 


87.8% 


82.6% 


50.0% 


9.2  xlO 


1.5  xlO' 


3.4  xlO' 


9.7  xlO' 


1  .4  xlO 


1.1  xlO 


1.54x0 


^II 

★ 

^1 

(Pii/Pg) 

10.3% 

10.6% 

32.0% 

32.9% 

30.0% 

31  . 9% 

23.8% 

2^7.1% 

18.76% 

22.7% 

12.5% 

16.6% 

1  . 76% 

3.5% 

DENSITY  FUNCTION 
OF  PREINSPECTION  FLAWS 
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same  phenomenon  is  observed  in  Table  7  that  the  Type  I  error  for  a  disk 
is  about  ten  orders  of  magnitude  smaller  than  that  for  a  hole  for 
mm.  Such  a  difference  will  be  manifested  drastically  as  the  number  of 
bolt  holes  increases.  A  comparison  between  Tables  6  and  7  indicates 
that  the  Weibull  POD  curve  is  more  beneficial  in  terms  of  the  Type  II 
error.  This  is  because  its  bandwidth  is  smaller.  However,  the  Type  I 
error  increases  more  rapidly  as  the  value  of  the  inspection  limit  aj^^^ 
reduces. 

c.  Exponential  POD  Function  with  a  Lower  Bound 

I 

;  During  the  first  few  inspection  maintenances,  the  majority  of 

flaws  are  still  very  small,  and  the  rejection  of  good  disks  (Type  II 
error)  is  of  main  concern,  as  demonstrated  in  the  previous  two  examples, 
j  Therefore,  we  shall  consider  a  POD  function  having  a  lower  bound  below 

which  the  crack  length  can  not  be  detected.  Let  the  POD  curve  be  of  the 
exponential  type  with  a  lower  bound  ag. 


P0D(a;l)  =  0 


a<ao 


=  1  -  exp{-6j^(a-aQ) }  a>aQ 


(34) 


in  which  a^  =  1  mm  and  6^=0.46/rmn.  Such  a  POD  curve  is  identical  to  the 
one  considered  in  the  first  example  except  that  the  curve  is  shifted  1 
mm  to  the  right-hand  side.  The  pre-inspection  flaw  distribution  is 
identical  to  that  of  the  previous  example. 

The  POD  curve  for  the  NOE  system  (Equation  34),  the  resulting 
POD  curve  for  a  disk,  and  the  pre-inspection  density  function  of  the 
flaw  length  are  displayed  in  Figure  18.  The  inspection  limit  correspond¬ 
ing  to  a  90ilJ  detection  probability  is  mm.  The  results  of  Type  I 

and  Type  II  errors,  the  conditional  Type  II  error,  and  the  average 
percentage  of  good  disks  are  presented  in  Table  8  for  different  values 
of  the  inspection  limit  a.inc*  ^  comparison  between  Tables  6  and  8 
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Figure  18.  Exponential  POD  Curves  with  Lower  Bound  at  1  mm  for  the 
NOE  System  and  Disk. 


TABLE  8 

TYPE  I  AND  TYPE  II  ERRORS:  EXPONENTIAL  POD  CURVE  WITH  LOWER 
BOUND  AT  a^  =  1  mm,  6^  =  0.46/mm 


No.  of 

Holes 

^NDE 

(mm) 

Pi 

Pii 

* 

PlI 

(Pii/Pq) 

1 

6 

98.5% 

9.1  xlO"^ 

18.3% 

18.5% 

10 

6 

98.5% 

WWIBH 

41.6% 

42.3% 

10 

5 

97.0% 

4.1  xlO'^^ 

40.0% 

41.2% 

10 

4 

94.0% 

37.9% 

39.3% 

10 

3 

nmi\m 

35.0% 

10 

2 

75 . 3% 

HSBByi 

18.4% 

24,4% 

10 

1  .5 

65.0% 

DSBBHI 

8.6% 

13.2% 
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indicates  that  the  POD  curve  with  a  lower  bound  does  improve  the  Type  II 
error  for  disks  whereas  the  Type  I  error  is  still  insignificant. 

d.  Narrow-Banded  Weibull  POD  Function 


For  the  purpose  of  illustration,  an  hypothetical  Weibull  POD 
curve  for  an  NDE  system  with  a  narrow  bandwidth  is  considered. 


in  which  aQ=8.0  and  6q=2.48  mm.  Such  a  narrow-banded  POD  curve  may  be 
achieved  by  an  NDE  system  when  the  bolt  hole  is  replicated.  The  POD 
curve  in  Equation  35  and  the  P0D(a;U10)  for  a  disk  are  displayed  in 
Figure  19.  The  pre-inspection  flaw  distribution  is  identical  to  that  of 
the  previous  example.  The  results  of  Type  I  and  Type  II  errors  and  the 
average  percentage  of  good  disks  are  presented  in  Table  9.  Again,  the 
Type  I  error  reduces  drastically  for  a  disk  as  compared  to  that  for  a 
hole.  The  Type  II  error,  however,  is  significantly  smaller  than  all  the 
previous  cases,  and  it  may  be  acceptable  for  the  RFC  system. 

e.  Conclusions 

The  conclusions  obtained  from  the  previous  examples  are  summarized 
below:  (1)  The  difference  in  the  Type  I  error  is  dramatic  between  inspect¬ 
ing  a  disk  with  ten  bolt  holes  and  inspecting  a  single  bolt  hole.  The 
difference  becomes  more  radical  as  the  number  of  bolt  holes  in  a  disk 
increases.  (2)  The  Type  I  error  for  disks  containing  many  bolt  holes 
appears  to  be  insignificant.  (3)  The  Type  II  error  for  disks  may  be  of 
practical  concern.  However,  if  the  bandwidth  of  the  POD  curve  of  the 
NDE  system  is  narrow,  as  indicated  by  the  last  example,  the  Type  II  error 
may  be  acceptable.  (4)  The  direction  for  improving  the  NDE  system  in 
applications  to  RFC  of  engine  components  appears  to  be  the  bandwidth  of 
the  POD  curve  rather  than  the  value  of  the  inspection  limit  that  cor¬ 
responds  to  a  90%  detection  probability  and  a  95%  confidence  level  as 
specified  by  damage  tolerance  requirements. 
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Figure  19.  Narrow-Banded  Weibull  POD  Curves  for  the  NDE  System 
and  Disk. 


TABLE  9 

TYPE  I  AND  TYPE  II  ERRORS:  NARROW-BANDED  WEIBULL  POD  CURVE  WITH 

Uq  =  8,0,  3q  =  2.48  mm 


No .  of 

Holes 

^NDE 

(mm) 

'’g 

'’l 

^i 

* 

Pii 

(Pii/Pq) 

1 

2.75 

85.4% 

4.4  xlO"^ 

3.8% 

4.4% 

10 

2.75 

85.4% 

15.3% 

17.9% 

10 

2.5 

82.6% 

BSBBI 

12.5% 

15.1% 

10 

2.25 

79.3% 

BWTWTHI 

9.2% 

11.0% 

10 

2.15 

77.8% 

BEE9SI 

7.7% 

9.9% 

10 

2.05 

76.2% 

IBHBI 

6.2% 

8.1% 

10 

2.0 

75.3% 

BEBBI 

7.3% 
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3.  TOTALLY  UNCORRELATED  CRACK  LENGTH 

At  another  extreme,  the  crack  length  varies  from  one  bolt  hole  to 
another  not  only  within  the  same  disk  but  also  across  the  ensemble  of 
disks.  The  crack  length  in  each  bolt  hole  is  a  statistically  independent 
and  identically  distributed  random  variable  as  shown  in  Figure  15(b). 

In  the  case  of  totally  correlated  crack  length,  the  average  percentage 
of  good  disks  is  identical  to  the  average  percentage  of  good  holes. 
However,  for  the  case  of  totally  uncorrelated  crack  length,  the  average 
percentage  of  good  disks  is  significantly  smaller  than  that  of  the  bolt 
holes.  To  appreciate  the  difference,  consider  100  bolt  holes,  each  with 
a  random  crack  length,  in  which  ten  of  them  are  bad  holes,  i.e.,  90% 
good  holes.  If  a  disk  contains  ten  bolt  holes  and  the  100  bolt  holes 
are  randomly  distributed  into  ten  disks,  then  the  number  of  good  disks 
ranges  from  nine  (all  the  bad  holes  are  put  into  one  disk)  to  zero  (one 
bad  hole  in  each  disk).  Therefore,  the  percentage  of  good  disks  ranges 
from  0%  to  90%.  As  a  result,  given  the  same  pre-inspection  distribution 
of  the  flaw  length,  the  return-to-service  interval  should  be  shorter  (or 
the  inspection  limit  aj^jQ^  should  be  larger)  to  maintain  the  same  level 
of  Type  I  error  as  the  previous  case. 

Let  a.  be  the  crack  length  in  the  jth  bolt  hole  of  a  disk.  Then 

sJ 

the  probability,  Pq,  that  a  disk  consisting  of  m  bolt  holes  is  a  good 
disk  is  given  by 

Since  a.  ( j=l,2, . . . ,m)  are  statistically  independent  and  identically 

J 

distributed  random  variables.  Equation  36  becomes 


P 


r: 


II  -  [F(a^„p)] 


(37) 
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in  which  F(a)  is  the  distribution  function  of  the  crack  length  in  a  hole 
prior  to  inspection. 

If  f(a)  denotes  the  probability  density  function  of  the  crack  length 
in  a  hole,  the  Type  I  error,  Pj,  for  a  disk  is  the  probability  of  accept¬ 
ing  a  bad  disk,  and  it  can  be  expressed  as 


{( 


f(x)[l-P0D(x;l )]dx 


NDE 


f(x)[l-P0D(x;l )]dx 


(38) 


in  which  the  first  integral  term  is  the  probability  of  accepting  a  bolt 
hole  (regardless  of  good  or  bad  hole),  and  hence,  the  first  term 
represents  the  probability  of  accepting  a  disk  or  m  bolted  holes  (again 
regardless  of  good  or  bad  disk).  The  second  integral  term  of  Equation 
38  is  the  probability  of  accepting  a  good  hole,  and  hence,  the  second 
term  indicates  the  probability  of  accepting  a  good  disk  or  m  good  holes. 

If  the  NDE  system  is  able  to  detect  all  the  cracks  larger  than  the 
inspection  limit  i.e.,  P00(a;l)=l  for  a2,a|^p£,  then  it  follows  from 

Equation  38  that  the  Type  I  error  is  zero  (i.e.,  Pj=0).  A  schematic 
flow  chart  presented  in  Figure  20  will  help  to  explain  Equations  38  and 
39. 

The  Type  II  error  for  a  disk  is  the  probability  of  rejecting  a  good 
disk,  which  can  be  expressed  as  follows; 


II 


=  Pr 


,  NDE 

U 


f(x)[l-POD(ax;l )]dx 


(39) 


in  which  Pg  =  [F(a|^[l^)]'^  given  by  Equation  37  is  the  probability  that 
the  disk  is  good.  The  second  term  on  the  right-hand  side  of  Equation  39 
is  the  probability  of  accepting  a  good  disk  (Figure  20).  If  the  NDE 
system  cannot  detect  any  crack  length  below  the  inspection  limit 
i.e.,  P0D(a;l)=0  for  a^ai^P^,  then  it  follows  from  Equation  39  that  the 
Type  II  error  is  zero,  i.e.,  Pjj=0. 
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Since  the  acceptance  of  a  disk  implies  the  acceptance  of  all  m  bolt 
holes  in  a  disk,  the  crack  length  of  a  disk,  denoted  by  a^,  can  be 
described  by  the  maximum  crack  length  of  m  holes,  i.e.. 


a 


m 


maxi 

l£j£m 


(40) 


Thus,  a  disk  is  good  if  its  crack  length  a^  is  smaller  than  the 
inspection  limit  The  distribution  function,  F^^^Ca),  of  the  crack 

length  a^  of  a  disk  can  be  derived  from  that  of  the  crack  length  of  a 
hole,  F(a),  using  the  extreme  value  distribution. 


F,„(a)  = 
m 


P[a  <a] 
m— 


■  m 

{a.<a} 

.j=l  ^  . 


=  [F(a)]" 


(41) 


The  probability  density  function  of  the  crack  length  a^^^  of  a  disk  is 
given  by  f^(a)=dFj^(a)/da  and  it  follows  from  Equation  41  that 

f^(a)  =  m[F(a)]"'*l  f(a)  (42) 


in  which  f(a)  is  the  probability  density  function  of  the  crack  length  in 
the  bolt  hole. 

Therefore,  the  probability  that  the  disk  is  good  is  given  by 


(431 


in  which  Equation  41  has  been  used.  As  expected,  Pg  derived  in  Equation 
43  is  identical  to  Equation  37. 

While  it  is  possible  to  describe  the  crack  length  of  a  disk,  a^, 
the  POD  curve  for  a  disk  loses  its  meaning  in  the  present  case  because 
of  the  following  reasons:  (1)  The  inspection  of  a  disk  implies  the 
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inspections  of  m  bolt  holes  in  which  the  crack  length  varies  independently 
from  one  hole  to  another;  and  (2)  The  inspection  of  a  disk  is  not 
identical  to  the  inspection  of  the  crack  length  a^^  of  the  disk.  As  a 
result,  the  POD  curve  for  a  disk  does  not  exist.  Consequently,  both  the 
capability  of  an  NDE  system  and  the  inspection  reliability  for  a  disk 
should  be  described  in  terms  of  Type  I  and  Type  II  errors. 

4.  EXAMPLE  No.  4 

a.  Various  POD  Functions 

For  the  purpose  of  comparison  between  the  two  extreme  cases 
investigated  previously,  the  same  numerical  example  for  the  case  of  totally 
correlated  crack  length  with  the  exponential  POD  curve,  and  the  distri¬ 
bution  of  the  pre-inspection  flaw  length  in  each  hole  displayed  in 
Figure  16,  are  considered.  Both  curves  are  also  shown  in  Figure  21  as 
the  dotted  curve  and  Curve  1,  respectively,  along  with  the  distribution 
of  the  pre-inspection  flaw  length  of  a  disk  (Equation  42)  shown  as  Curve 
2.  It  can  be  observed  from  Figure  21  that  the  distribution  of  the  pre¬ 
inspection  flaw  length  for  a  disk  is  shifted  to  the  right-hand  side, 
indicating  a  larger  percentage  of  bad  disks  than  that  of  bad  holes. 


Figure  21.  Distributions  of  Pre-inspection  Flaw  Length  for  Bolt 
Hole  and  Disk;  x  =  0.7/mm. 
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The  average  percentage  of  good  disks,  Pg,  is  obtained  by  sub¬ 
stituting  Equation  14  into  Equation  43, 

Pg  =  [1  -  exp  (< 

and  the  results  are  given  in  Table  10  for  different  values  of  the 


inspection  limit  a. 


A  comparison  between  Tables  6  and  10  indicates  a 


substantial  reduction  of  the  average  percentage  of  good  disks  prior  to 
inspection,  when  the  crack  length  in  each  hole  is  statistically 
independent. 


TABLE  10 

TYPE  I  AND  TYPE  II  ERRORS  FOR  TOTALLY  UNCORRELATED  CRACK  LENGTH; 
EXPONENTIAL  POD  CURVE  WITH  3,  =  0.46/mm  and  X  =  0.7/mm 


No.  of 
Holes 


Pi 

PlI 

★ 

''ii 

(Pii/pg) 

97.0% 

1. 830x10"^ 

36.8% 

37.9% 

73.6% 

BUfTHTM 

72.9% 

99.0% 

53.4% 

52.8% 

98.9% 

27.1% 

26.6% 

98.2% 

86.0% 

85.3% 

99.2% 

92.8% 

92.2% 

99.4% 

96.4% 

95.72% 

99.3% 

98.2% 

0.187x10'^ 

97.5% 

99.29% 
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The  Type  I  and  Type  II  errors  can  be  obtained  by  substituting 
Equations  14  and  30  into  Equations  38  and  39;  with  the  results 


in  which  Pg  is  given  in  Equation  44. 

The  results  of  Type  I  and  Type  II  errors  as  well  as  the  con¬ 
ditional  Type  II  error  are  presented  in  Table  10  for  different  values  of 
the  inspection  limit 


Although  the  Type  I  error  at  5  mm  for  a  disk  is  smaller 

than  that  for  a  hole  in  Table  10,  the  improvement  is  by  no  means  as  sig¬ 
nificant  as  that  shown  in  Table  6.  It  is  impossible  to  achieve  the  kind 
of  Type  I  error  shown  in  Table  6,  even  when  the  inspection  limit  aj^p^  is 
increased  to  9  mm.  While  the  Type  II  error,  Pjj,  for  a  disk  at  aj^p^  =  5  mm 
is  72.9%,  the  conditional  Type  II  error,  Pj^,  is  99%  because  the  average^ 
percentage  of  good  disks  is  only  73.6%.  In  other  words,  if  we  know  that 
the  disk  is  good,  the  probability  of  rejecting  it  is  99%.  The  last 
column  of  Table  10  indicates  that  the  Type  II  error  is  extremely  bad 
when  the  POD  curve  of  the  NDE  system  is  the  exponential  function  and  the 
pre- inspect  ion  crack  length  in  each  hole  is  statistically  independent. 


In  the  first  few  inspection  maintenances,  the  average  percentage 
of  bad  disks  is  usually  very  small  because  of  safety  considerations. 

Thus,  the  example  considered  previously  may  be  more  representative  when 
the  disk  is  in  service  for  a  longer  period  of  time.  We  shall  investigate 
the  situation  in  which  the  pre-inspection  flaw  length  at  holes  is  reduced. 
In  this  connection,  the  same  exponential  distribution.  Equation  14,  will 
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be  used,  but  with  a  smaller  average  flaw  length,  i.e.,  l/x  =  1.0  mm 
or  A  =  1.0/mm,  as  shown  in  Figure  22  by  Curve  1.  Also  plotted  in  Figure 
22  as  Curve  2  is  the  distribution  of  the  pre-inspection  flaw  length  for 
a  disk.  The  average  percentage  of  good  disks  for  a|yjQ^=  5  mm  prior  to 
inspection  is  computed  from  Equation  44  as  93.46%.  The  same  POD  curve 
employed  in  the  previous  example  is  used  herein  and  plotted  as  a  dashed 
curve  in  Figure  22. 

The  results  of  the  average  percentage  of  good  disks.  Type  I  and 
Type  11  errors,  and  conditional  Type  II  error  are  presented  in  the  upper 
part  of  Table  11  for  different  values  of  the  inspection  limit 
Unlike  Table  6,  the  Type  I  error  in  Table  11  improves  only  slightly  at 
^NDE“  ^  ^  deteriorates  as  the  inspection  limit 

reduces.  Likewise,  both  the  Type  II  error  and  the  conditional  Type  II 
error  are  too  high  for  practical  applications. 

Instead  of  the  exponential  POD  curve  employed  above,  the  Weibull 
type  POD  curve  shown  in  Figure  17  as  P00{a;l),  i.e.,  oq  =  3.0  and  = 

3.79  mm,  is  considered.  The  results  of  Type  I  error.  Type  II  error,  and 
conditional  Type  II  error  are  presented  in  the  lower  part  of  Table  11. 

It  is  observed  that  the  Type  I  error  for  a  disk  is  worse  than  that  for  a 
hole  at  5  mm,  and  it  becomes  much  worse  as  the  inspection  limit 

^NDE  However,  the  Type  II  error  improves  over  that  of  the 

previous  example,  although  it  is  still  quite  significant. 

Finally,  let  us  consider  a  narrow-banded  Weibul 1-type  POD  curve 
as  shown  by  Curve  1  of  Figure  12,  i.e.,  Og  =  8.0  and  6q  =  4.13  mm.  The 
results  are  shown  in  Table  12.  Table  12  demonstrates  a  significant 
improvement  over  Table  11,  indicating  the  importance  of  the  bandwidth  of 
the  POD  curve.  In  Table  12,  both  Type  I  and  Type  II  errors  may  be 
acceptable.  In  reality,  the  narrow-band  POD  curve  used  in  this  example, 
i.e.,  oiq  =  8.0  and  3q  =  4.13  mm,  may  not  be  achieved  by  an  NDE  system. 
Fortunately,  however,  such  a  narrow-banded  POD  curve  can  be  established 
by  two  inspections  using  the  intersection  rule  described  previously  in 
which  the  rejected  holes  may  be  cleaned  up,  polished  or  even  replicated 
before  the  second  inspection,  as  described  in  the  next  section. 
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CRACK  LENGTH,  MM 

Figure  22.  Distribution  of  Pre-inspection  Flaw  Length  for  Bolt 
Hole  and  Disk;  A  =  1.0/mm. 


TABLE  11 

TYPE  I  AND  TYPE  II  ERRORS  FOR  TOTALLY  UNCORPELATED  CRACK  LENGTH 

WITH  A  =  1 . O/mm 

EXPONENTIAL  TYPE  POD  CURVE;  3,  =  0. 46/mm 


No.  of  Pg 

Holes  (mm) 


5  99.33%  4.60x10" 


93.46%  I  1 .53x10 


83.12%  I  6.50x10 


60.00%  2.60x10" 


WEIBULL  TYPE  POD  CURVE;  a.  =  3.C 


99.33%  I  2.60x10 


93.46%  1.40x10" 


83.12%  1.42x10 


60.00%  8.90x10" 


0.46/mm 

'll 

* 

'll 

(Pll/Pg) 

30.9% 

Hffli 

91 .2% 

97.6% 

80.9% 

97.3% 

58.0% 

97.0% 

o 

11 

CO 

79  mm 

6.1% 

6.14% 

43.7% 

46.76% 

3^ .  6% 

41.6% 

18.9% 

31  .5% 
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TABLE  12 


TYPE  I  AND  TYPE  II  ERROR  FOR  TOTALLY  UNCORRELATED  CRACK  LENGTH 
WITH  X  =  1.0/mm  AND  NARROW-BANDED  WEIBULL-TYPE  POD  CURVE 
FOR  aQ  =  8.0  AND  3q  =  4.13  mm 


No.  of 

Holes 

^NDE 

(mm) 

Pl 

PlI 

'll 

(Pll/Pg) 

1 

5 

99.332 

7.00x10“^ 

1  . 79% 

1  .8% 

10 

5 

93.46% 

BB8ISI 

15.50% 

16.6% 

10 

4.8 

92.07% 

14.20% 

15.4% 

10 

4.6 

90.39% 

12.6% 

13.9% 

10 

88.38% 

3.71x10"^ 

10.8% 

12.2% 

b.  Conclusions 

The  conclusions  obtained  previously  for  the  totally  uncorrelated 
crack  length  are  summarized  as  follows:  (1)  The  average  number  of  bad 
disks  is  much  higher  than  the  average  number  of  bad  holes  prior  to  in¬ 
spection.  This  is  of  practical  importance  because  in  the  Retirement- 
For-Cause  analysis,  the  distribution  of  the  flaw  length  at  each  hole  is 
simulated.  (2)  Both  Type  I  and  Type  II  errors  are  more  serious  than 
those  in  the  case  of  totally  correlated  crack  length.  (3)  Depending  on 
the  pre-inspection  flaw  length  in  holes  or  equivalently  the  length  of 
service  life  for  disks  prior  to  inspection,  both  Type  I  and  Type  II  errors 
may  be  of  concern,  if  the  bandwidth  of  the  POD  curve  is  not  narrow  enough. 
(4)  When  the  POD  curve  is  narrow-banded,  both  Type  I  and  Type  II  errors 
may  be  acceptable,  as  demonstrated  in  the  last  numerical  example. 
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SECTION  V 

MULTIPLE  INSPECTIONS  FOR  ENGINE  DISKS 

It  was  explained  in  the  previous  section,  that  if  the  bandwidth  of 
the  POD  curve  of  the  NDE  system  is  not  narrow  enough,  then  the  Type  II 
error  may  be  of  great  concern  in  the  RFC  system.  Should  this  be  the 

case,  a  redundant  NDE  system  for  a  second  inspection  of  rejected  disks 

is  highly  desirable,  as  will  be  demonstrated  in  this  section.  Mathe¬ 
matical  results  for  engine  disks  are  presented,  in  which  each  hole 

rejected  by  the  No.  1  NDE  system  is  further  inspected  by  the  No.  2  NDE 

system  (intersection  rule).  Numerical  examples  are  given  to  illustrate 
quantitatively  the  advantage  of  multiple  inspections  for  engine  disks. 

To  effectively  reduce  the  Type  II  error,  it  is  assumed  that  the  POD 
curve  of  the  No.  2  (or  redundant)  NDE  system  is  narrow-banded,  represent¬ 
ing  the  inspection  of  holes  with  polishment  or  replication.  Again,  two 
extreme  cases  will  be  investigated;  totally  correlated  flaw  length  and 
totally  uncorrelated  flaw  length. 

1.  TOTALLY  CORRELATED  FLAW  LENGTH 

When  a  hole  is  inspected  by  two  NDE  systems  with  P0D(a;l)  and 
P0D(a;2),  respectively,  using  the  intersection  rule,  the  resulting  POD 
curve  for  a  hole  follows  from  Equation  4  as 

P0D(a;la2)  =  P0D(a;l)P0D(a;2)  (47) 

The  POD  curve  for  a  disk  consisting  of  m  holes,  denoted  by  P0D(a;Um), 
follows  from  Equations  27  and  47  as 

P0D(a;Um)  =  1  -  [l-P0D(a;l)P0D(a;2)]"’  (48) 

Hence,  the  Type  I  and  Type  II  errors  for  disks  is  given  by 

'D 

Pj  =  J  f(x)[l-P0D(x;l  )P0D(x;2)f  dx  (49) 

'’nde 
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and 

^NDE 

Pj!  =  J  f(x)|l  -  [l-POD(x;l  )P0D(x;2)f|dx  (50) 

0 

in  which  f(a)  is  the  probability  density  function  of  the  pre-inspection 
flaw  length  in  holes.  The  probability,  Pg,  that  a  disk  is  good  (or  the 
average  percentage  of  good  disks)  is  given  by  Equation  28. 

2.  TOTALLY  UNCORRELATEO  FLAW  LENGTH 

Both  Type  I  and  Type  II  errors  as  well  as  the  average  percentage  of 
good  disks  are  derived  in  Equations  38,  39  and  37,  respectively,  in  which 
P0D(x;l)  should  be  replaced  by  P0D(x;l)P0D(x;2) ,  i.e., 

Pg  '  [F(aNDE>j"'  <5" 


A  flow  chart  explaining  the  procedures  of  two  inspections  using  the 
intersection  rule  for  disks  derived  above  is  given  in  Figure  23.  It 
should  be  mentioned  that  only  those  holes  rejected  by  the  No.  1  NDE  system 
are  further  inspected  by  the  No.  2  NDE  system  (Figure  23).  Therefore, 
the  probability,  Pn  ,  that  a  hole  will  be  rejected  by  the  No.  1  NDE 
system  is  given  by 
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ONEDISK 

■ 

NO.  1  NDE  SYSTEM 

m-m^  HOLES 

m  HOLES 

POD  (a;  1) 

ACCEPTED 

Figure  23.  Procedures  for  Two  Inspections  for  Disks. 
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Furthermore, 
inspected  by 


pQ  also  represents  the  average  percentage  of  holes  to  be 

the  No.  2  NDE  system.  Indeed,  Pn  is  an  important  quantity 

*^1 


in  determining  the  cost  benefit  for  two  inspections.  It  is  interesting 
to  observe  from  Equation  54  that  the  average  percentage  of  holes  to  be 
inspected  by  the  No.  2  NDE  system  is  independent  of  the  inspection  limit 

®NDE' 

3.  NUMERICAL  EXAMPLES 


As  mentioned  previously,  the  redundant  system  has  a  high  resolution 
capability,  representing  the  inspection  of  holes  with  polishment  or 
replication.  Therefore,  a  narrow-banded  Weibull  POD  curve  is  considered 


P0D(a;2)  =  1  -  exp[-(a/e)“]  ;  a>0  (55) 


in  which  a=8.0  and  e=4.13  mm.  Such  a  POD  curve  was  shown  in  Figure  12 
as  Curve  1. 

The  exponential  distribution  for  the  pre-inspection  flaw  length  given 
by  Equation  14  is  considered, 

f(a)  =  xe"^^  ;  a^O  (56) 


a.  Examples  for  Totally  Correlated  Flaw  Length 

For  the  case  of  totally  correlated  flaw  length,  we  choose  X=0.7/mm 
in  Equation  56.  Such  a  distribution  of  the  pre-inspection  flaw  length 
was  shown  in  Figure  16.  For  the  purpose  of  comparison  between  the  results 
of  single  inspection  and  that  of  two  inspections,  the  same  P0D(a;l)  shown 
in  Figure  16  is  considered,  i.e.,  P0D(a;l)  =  1  -  exp(-Bj^a)  with  62^=0. 46/mm. 

The  results  of  Type  I  error.  Type  II  error  and  conditional  Type 
II  error  for  disks  are  presented  in  Table  13.  The  corresponding  results 
under  a  single  inspection  were  shown  in  Table  6.  A  comparison  between 
Tables  6  and  13  indicates  a  drastic  reduction  for  the  Type  II  error 
resulting  from  the  second  inspection,  whereas  the  Type  I  error  remains 
insignificant.  For  instance,  with  the  inspection  limit 
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TABLE  13 


TYPE  I  AND  TYPE  II  ERRORS  FOR  TOTALLY  CORRELATED  CRACK 
LENGTH  USING  TWO  INSPECTIONS;  X  =  0.7/mtn,  6,  =  0.46/mm, 
a  =  8.0,  3  =  4.13  mm 


No.  of 

Holes 

^NDE 

(mm) 

’’g 

”1 

PlI 

★ 

PlI 

★ 

1 

5 

97.00% 

1 .80x10"^ 

36.8% 

37.9% 

1 

5 

97.00% 

1.90x10"^ 

3.5% 

3.6% 

10 

5 

97.00% 

6.04x10'^^ 

9.6% 

9.9% 

10 

4.5 

95.70% 

1 .22x10"^ 

8.3% 

8.7% 

10 

4.0 

94.00% 

7.45x10"^ 

6.5% 

6.9% 

10 

3.5 

91.37% 

1. 10x10"^ 

4.1% 

4.5% 

10 

3.0 

87.80% 

1 .36x10"^ 

1.7% 

1.9% 

* 

One  Inspection  For  One  Hole 

Type  II  error  under  a  single  inspection  is  80.7%,  whereas  it  is  reduced 
to  6.5%  using  two  inspections.  Likewise,  the  Type  I  errors  for  both 
cases  are  acceptable. 

The  average  percentage  of  holes,  P^  ,  to  be  inspected  by  the 

No.  2  NDE  systems  is  computed  from  Equation  54  as  40%.  Thus,  by  inspecting 
40%  of  the  holes  using  the  No.  2  NDE  system,  one  can  reduce  the  Type  II 
error  for  disks  from  80.7%  to  6.5%. 

As  a  second  example,  the  Weibull-type  P0D(a;l)  function  given 
by  Equation  33  with  aQ=3.0  and  3q=3.79  mm  is  considered.  Such  a  POD 
curve  was  shown  in  Figure  17.  The  results  of  Type  I  error.  Type  II  error 
and  conditional  Type  II  error  for  disks  are  shown  in  Table  14. 

It  is  observed  from  Table  14  that  the  Type  I  error  is  insig¬ 
nificant  and  the  Type  II  error  is  acceptable.  The  corresponding  results 
under  a  single  inspection  were  given  in  Table  7.  A  comparison  between 
Tables  7  and  14  indicates  that  a  significant  reduction  of  Type  II  error 
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can  be  achieved  by  the  second  inspection.  The  average  percentage  of 
holes  to  undergo  the  second  inspection  is  computed  from  Equation  54  as 
13%.  For  instance,  with  the  inspection  limit  mm,  the  Type  II 

error  is  reduced  from  30%  (Table  7}  to  4.7%  (Table  14)  by  inspecting  13% 
of  holes  (on  the  average)  for  a  second  time. 

TABLE  14 

TYPE  I  AND  TYPE  II  ERRORS  FOR  TOTALLY  CORRELATED  CRACK  LENGTH  USING  TWO 
INSPECTIONS;  X  =  0.7/mm,  aQ  =  3.0,  Bq  =  3.79  mm, 

a  =  8.0.  6  =  4,13  mm 


No .  of 


2.96% 


87.80%  2.16x10 


8.06% 


6.84% 


4.76%  5.06% 


2.47%  I  2.70% 


0.70%  0.80% 


One  Inspection  For  One  Hole. 


b.  Examples  for  Totally  Uncorrelated  Flaw  Length 

For  the  case  of  totally  uncorrelated  flaw  length,  x=1.0/mm  is 
chosen  for  Equation  56  to  compare  the  results  with  those  under  a  single 
inspection.  The  distribution  of  the  pre-inspection  flaw  length  was  shown 
in  Figure  22  as  Curve  1.  The  same  POD  curves  used  in  two  previous  examples 
i.e.,  P0D(a;l)  =  1  -  exp(-Bja)  with  Bj=0. 46/mm  and  P0D(a;l)  =  1  -  exp[- 
(a/BQ)“0]  (Equation  33)  with  0^=3. 0  and  6q=3. 79/mm  are  considered.  The 
results  for  Type  I  error.  Type  II  error  and  conditional  Type  II  error 
are  shown  in  Table  15.  The  corresponding  results  using  a  single  inspection 
were  given  in  Table  11.  A  comparison  between  Tables  11  and  15 


AFWAL-TR-82-4111 


demonstrates  a  significant  reduction  for  the  Type  II  error  using  a  second 
inspection  (Equation  55). 

TABLE  15 


4.  CONCLUDING  REMARKS 

When  the  bandwidth  of  the  POD  curve  is  not  narrow  enough,  the  second 
inspection,  using  an  NOE  system  having  a  narrow-banded  POD  curve,  will 
reduce  the  Type  II  error  substantially.  This  conclusion  holds  for  both 
cases  when  the  flaw  length  in  each  hole  of  the  same  disk  is  completely 
correlated  or  completely  independent.  The  second  inspection  may  be 
expensive  and  time  consuming,  for  instance,  the  replication  of  holes. 
However,  the  average  percentage  of  holes  to  undergo  the  second  inspection 
is  reasonably  small  and  therefore  the  second  inspection  may  be  economically 
beneficial . 
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Although  the  flaw  length  in  disk  holes  and  that  for  the  POD  curves 
used  in  the  numerical  examples  are  not  close  to  the  situation  for  gas 
turbine  engine  disks,  the  conclusions  and  trends  obtained  for  both  Type 
I  and  Type  II  errors  remain  the  same.  This  is  because  both  Type  I  and 
Type  II  errors,  as  well  as  the  average  percentage  of  good  holes  prior  to 
inspection,  remain  identical  when  the  flaw  length  for  POD  curves  and 
that  in  disk  holes  are  reduced  by  a  factor  of  y  (e.g.,  y  =  5.0). 
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SECTION  VI 


INSPECTION  RELIABILITY  OF  ENGINE  DISKS  USING  P&W  POD  CURVE 


The  numerical  examples  presented  in  the  previous  two  sections  are 
intended  only  to  illustrate  the  effect  of  various  types  of  POD  curves  on 
the  inspection  reliability  as  well  as  the  effect  of  multiple  inspections. 
While  the  POD  curves  used  are  obtained  from  the  literature,  they  are  far 
from  current  engine  inspection  technology,  although  the  trends  and  con¬ 
clusions  regarding  both  types  of  errors  will  not  be  altered.  Recently,  a 
laboratory  POD  curve  was  established  by  the  Pratt  and  Whitney  Aircraft 
Company,  as  shown  in  Figure  24.  Such  a  POD  curve  may  be  representative 
of  future  inspection  technology  for  RFC  inspection  systems,  and  it  is 
used  in  this  section.  Therefore,  examples  are  given  to  demonstrate  the 
magnitude  of  Type  I  and  Type  II  errors  for  engine  disks  under  a  single 
inspection,  but  using  different  distribution  functions  for  the  pre¬ 
inspection  flaw  length.  The  application  of  multiple  inspections  is  not 
given  due  to  lack  of  a  POD  curve  for  holes  with  polishment  or  replication. 

The  POD  curve  shown  in  Figure  24  can  be  fitted  by  the  three-parameter 
Wei  bull  distribution  very  well. 


P0D(a;l )  =  0 


for  a_<aQ 
for  a>aQ 


(57) 


in  which  ag  is  the  lower  bound  below  which  the  crack  length  can  not  be 
detected.  The  best-fitted  parameter  values  are  obtained  as  follows: 

_3 

ag  =  1,4  mils  (10  inch),  6^  =  '^•^5  mils  and  =  2.0.  With  this  set 
of  parameter  values,  the  plot  of  Equation  57  cannot  be  distinguished 
from  the  POD  curve  shown  in  Figure  24. 

The  negative  exponential  distribution  for  the  pre-inspection  flaw 
length  given  by  Equation  14  is  employed.  Currently,  the  target  for  the 

_  3 

inspection  limit,  of  advanced  engine  dirks  is  about  5  mils  (5x10 

inches).  Consequently,  the  average  pre-inspection  flaw  length  is  assumed 
to  be  1/X  =  1.429  mils,  or  X  =  0.7/mil,  indicating  that  at  the  time  of 
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inspection,  3%  of  holes  will  have  a  flaw  length  longer  than  5  mils,  i.e., 
Pg  =  0.97.  The  density  function  f(a)  of  the  pre-inspection  flaw  length 
is  shown  in  Figure  25  as  Curve  1. 


The  results  of  Type  I  and  Type  II  errors  are  presented  in  Table  16, 
in  which  "+"  represents  the  case  for  totally  correlated  flaw  length  in 
each  hole  of  the  same  disk  and  "++"  denotes  the  case  for  totally 
independent  flaw  length  in  each  hole.  Table  16  indicates  similar  trends 
already  shown  in  Tables  6-12,  as  follows. 


(1)  For  the  case  of  completely  correlated  flaw  length,  the  Type  I 
error  is  very  small.  The  Type  II  error,  however,  is  moderate  (e.g.,  at 
^NDE  ~  ^  HI  is,  Pjj  =  15.3%).  This  is  so  because  the  POD  curve  has  a 
lower  bound  at  ag  =  1.4  mils  where  62.5%  of  the  flaws  are  shorter  than 
1.4  mils.  The  beneficial  effect  resulting  from  a  POD  curve  which  has  a 
lower  bound  has  been  discussed  previously. 


(2)  For  the  case  of  totally  independent  flaw  length,  the  Type  I 
error  requires  a  serious  consideration  and  the  Type  II  error  is  also 
higher  in  the  range  where  mils.  For  instance,  at  a|^g^=5  mils,  we 


NDE- 

have  Pjj  =  23.4%  and  Pjj  =  31.8%. 


Suppose  the  return  to  service  interval  is  increased,  or  the  safety 
factor  is  reduced,  such  that  at  the  time  of  inspection,  there  are  8.2%  of 
flaws  exceeding  5  mils,  i.e.,  x  =  0.5/mi  1  or  the  average  flaw  length 
=  1/A  =  2  mils.  Such  a  distribution  is  shown  in  Figure  25  as  Curve  2. 

The  results  of  Type  I  and  Type  II  errors  are  presented  in  Table  17.  For 
the  case  of  totally  correlated  flaw  length,  the  Type  I  error  is  always 
small  but  the  Type  II  error  increases  over  the  previous  example,  as 
expected,  because  of  more  large  flaws.  For  the  case  of  totally  independent 
flaw  length,  the  Type  I  error  requires  considerable  attention.  Although 
the  Type  II  error  at  a|^g^  =  5  mils  is  Pjj  =  19.4%,  the  conditional  Type  II 
error  increases  to  Pjj  =  45.7%.  This  is  because  we  have  only  42.5%  of 
good  disks  at  the  time  of  inspection,  i.e.,  Pp  =  42.5%. 
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TABLE  16 


TYPE  I  AND  TYPE  II  ERRORS  FOR  DISK  USING  P&W  POD  CURVE  AND 
AVERAGE  FLAW  LENGTH  1,429  mils  (X=0. 71 /mi  1 ) ;  EXPONENTIAL 
DISTRIBUTION  FOR  PRE-INSPECTION  FLAW  LENGTH,  +  TOTALLY 
CORRELATED  FLAW  LENGTH,  ++  TOTALLY  INDEPENDENT  FLAW  LENGTH 


NO.  OF 

HOLES 

^G 

1 

93.9% 

lo" 

93.9% 

10^" 

53.4% 

1 

97.0% 

lO"^ 

97.0% 

lo"^"^ 

73.6% 

1 

98.5% 

lo" 

98.5% 

10^^ 

86.0% 

1 

99.3% 

10^ 

99.3% 

lO'^'^ 

92.8% 

1 

99.6% 

lo" 

99.6% 

10^^ 

96.4% 

1 

99.8% 

10^ 

99.8% 

lo"^"^ 

98.2% 

7.79x10' 

1.47x10' 

4.35x10' 


2.15x10' 

4.04x10' 

1.23x10' 


1.15x10' 

8.53x10' 

6.67x10' 


2.24x10' 

6.49x10' 

1.29x10' 


2.3% 

12.3% 

11.8% 


3.6% 

15.3% 

23.4% 


4.6% 

16.9% 

32.7% 


5.2% 

17.6% 

38.6% 


5.5% 

17.9% 

41.9% 


5.7% 

18.2% 

43.7% 


2.5% 

13.1% 

22.1% 


3.7% 

15.8% 

31.8% 


4.7% 

17.2% 

38.0% 


5.2% 
17.7% 
41 .6% 


5.5% 

18.0% 

43.5% 


5.7% 

18.2% 

44.5% 
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TABLE  17 

TYPE  I  AND  TYPE  II  ERRORS  FOR  DISK  USING  P&W  POD  CURVE  AND 
AVERAGE  FLAW  LENGTH  2.0  mils  (A=0.5  mil);  EXPONENTIAL 
DISTRIBUTION  FOR  PRE-INSPECTION  FL^W  LENGTH,  +  TOTALLY 
CORRELATED  FLAW  LENGTH,  ++  TOTALLY  INDEPENDENT  FLAW  LENGTH 


B 

NO.  OF 

HOLES 

Pg 

Pl 

PlI 

★ 

^II 

■■ 

1 

86.5% 

4.79x10'^ 

3.1% 

3.6% 

4 

lO"^ 

86.5% 

2.77x10'^ 

15.8% 

18.3% 

10^" 

23.4% 

1.21x10"’ 

7.2% 

20.8% 

1 

CO 

1.78x10"^ 

5.4% 

5.9% 

5 

lO"^ 

91 .8% 

2.97x10"^ 

21.2% 

23.0% 

10-^" 

42.5% 

5.23x10"^ 

19.4% 

45.7% 

1 

95.0% 

5.90x10"^ 

7.5% 

7.9% 

6 

lo"^ 

95.0% 

9.87x10"^ 

24.4% 

25.7% 

10^" 

60.0% 

1.84x10"^ 

33.5% 

55.8% 

1 

97.0% 

1.75x10"^ 

9.0% 

9.3% 

7 

lo" 

97.0% 

1.00x10"” 

26.3% 

27.1% 

73.6% 

5.58x10"^ 

45.8% 

62.2% 

1 

98.2% 

4.63x10"^ 

10.1% 

10.3% 

8 

lo"^ 

98.2% 

3.09x10"’^ 

27.5% 

28.0% 

83.1% 

1 .49x10"^ 

55.0% 

66.2% 

1 

98.9% 

1 .09x10"^^ 

10.8% 

10.9% 

9 

lo"^ 

98.9% 

2.86x10"’^ 

28.2% 

28.5% 

10^" 

89.4% 

3.50x10"^ 

61.2% 

68.5% 
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For  all  cases  considered  so  far,  including  those  of  the  previous 

★ 

sections,  the  conditional  Type  II  error  for  the  case  of  completely 
independent  flaw  length  is  always  higher  than  that  for  the  case  of 
totally  correlated  flaw  length. 

As  mentioned  previously,  the  negative  exponential  distribution  for 
the  pre-inspection  flaw  length  given  by  Equation  14  is  a  special  case  of 
the  Weibull  distribution  in  which  the  coefficient  of  variation  (dispersion 
is  100/!;.  It  may  be  interesting  to  examine  the  sensitivity  of  Type  I  and 
Type  II  errors  with  respect  to  the  variation  of  the  distribution  function 
used  for  the  pre-inspection  flaw  length.  In  this  connection,  the  log¬ 
normal  distribution  given  by  Equation  14(b),  is  considered.  The 
coefficient  of  variation,  V,  of  the  lognormal  distribution  is  related 
to  a  through 

V  =  ■|exp[(ailn  10)  ]  -  1  (58) 

To  be  consistent  with  the  previous  examples,  V  is  chosen  to  be  100%, 
resulting  in  a  =  0.3616.  Two  cases  corresponding  to  the  previous  examples 
for  A  =  0.7/mil  and  A  =  0.5/mil,  respectively,  are  considered:  (1)  a 
median  flaw  length  of  1.045  mils  (i.e.,  u  =  0.019),  that  results  in  3%  of 
flaws  exceeding  a^^^^  =  5  mils,  and  (2)  a  median  flaw  length  of  1.57  mils 
(i.e.,  u  =  0.1963),  that  results  in  8.2%  flaws  exceeding  a^^^^  =  5  mils. 
These  two  lognormal  probability  density  functions  are  shown  in  Figure  25 
as  Curve  3  (u  =  0.019)  and  Curve  4  (p  =  0.1963),  respectively,  for 
compari son . 

The  results  of  Type  I  and  Type  II  errors  are  presented  in  Tables  18 
and  19,  respectively,  for  u  =  0.019  and  y  =  0.1963.  A  comparison  between 
Tables  16  and  18  as  well  as  Tables  17  and  19  indicates  that  the  difference 
is  very  small.  This  is  expected  from  Figure  24,  since  the  main  difference 
between  the  exponential  and  lognormal  distributions  lies  in  the  small 
flaw  size  region,  a<ag  =  1.4  mil,  in  which  the  NOE  system  cannot  detect 
any  flaw.  Consequently,  when  the  number  of  bad  holes  is  within  8%  at  the 
time  of  inspection,  the  functional  form  for  the  distribution  function  of 
the  pre-inspection  flaw  length  is  not  important,  as  long  as  the 
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■,t.jtistical  dispersion  is  identical  and  there  is  a  reasonable  lower 
bound  for  the  POD  curve. 


TABLE  18 

TYPE  I  AND  TYPE  II  ERRORS  FOR  DISK  USING  P&W  POD  CURVE  AND 
MEDIAN  FLAW  LENGTH  1.045  mils  (p=0.019);  LOGNORMAL  DISTRI¬ 
BUTION  FOR  PRE- INSPECTION  FLAW  LENGTH,  +  TOTALLY  CORRELATED 
FLAW  LENGTH,  ++  TOTALLY  INDEPENDENT  FLAW  LENGTH. 


H 

NO.  OF 

HOLES 

^G 

>^1 

'’ll 

•k 

1 

94.7% 

2.0% 

2.1% 

4 

10"^ 

94.7% 

10.7% 

11.3% 

10^^ 

57.7% 

1 .10x10'' 

10.9% 

18.9% 

1 

97.0% 

6.50x10'^ 

3.0% 

5 

lo"^ 

97.0% 

1.15x10'^ 

13.0% 

IC^^ 

73.7% 

3.86x10'^ 

19.7% 

26.7% 

1 

98.2% 

2.01x10'^ 

3.7% 

3.8% 

b 

lO"" 

98.2% 

3.45x10'^ 

14.2% 

14.5% 

10^" 

83.5% 

1.22x10'^ 

26 . 8% 

32.1% 

1 

98.9% 

5.82x10'^ 

4.4% 

7 

10^ 

98.9% 

3.33x10"'^ 

BH 

15.1% 

10^" 

89 . 3% 

3.55x10'^ 

31  . 7% 

35.5% 

1 

99.3% 

1 .56x10'^ 

4.6% 

4.6% 

8 

lo”^ 

99.3% 

1 .02x10"'^ 

15.3% 

15.4% 

lo"-^ 

93.0% 

9.53x10'^ 

35 . 2% 

37.9% 

1 

99.5% 

3.84x10'^ 

4.8% 

4.8% 

-20 

9 

10 

99.5% 

9.72x10 

15.5% 

15.6% 

10^^ 

95.2% 

2.35x10'^ 

37.4% 

39.3% 
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Itr.l.f'.nl  Ilf  ■l^^llllli  ri')  .)  1007.  coefficient  of  variation  (dispersion) 

for  the  pre-inspection  flaw  length,  a  507  coefficient  of  variation  is 

considered  in  conjunction  with  the  lognormal  distribution  given  by 

Equation  14(b).  The  parameter  a  is  then  computed  from  Equation  58  with 

V  =  507  as  0.2052,  i.e.,  a  =  0.2052.  For  the  purpose  of  comparison  with 

the  previous  example,  the  same  percentage  of  flaws  smaller  than 

mils  at  the  time  of  inspection  is  considered,  i.e.,  Pg  =  977.  Then  u  is 

calculated  from  Eq.  14(b)  as  0.3133,  indicating  that  the  median  flaw 
0  31  33 

length  is  (10)  '  =  2.06  mils.  Such  a  lognormal  probability  density 

function  for  the  pre-inspection  flaw  length  is  displayed  in  Figure  26. 

A  comparison  between  Curve  3  of  Figure  25  and  Figure  26  shows  that 
when  the  statistical  dispersion  of  the  pre-inspection  flaw  length  is 
reduced,  the  peak  of  the  density  function  is  shifted  to  the  right-hand 
side  (large  crack  size  region)  in  order  to  maintain  the  same  percentage 
of  good  holes.  As  a  result,  it  is  expected  that  the  Type  II  error  will 
i ncrease. 

The  Type  I  and  Type  II  errors,  the  percentage  of  good  disks,  and  the 
conditional  Type  II  error  are  shown  in  Table  20.  A  comparison  between 
Tables  18  and  20  indicates  that  the  Type  II  error  increases  substantially 
for  the  inspection  limit  mils.  For  instance,  with  mils, 

the  Type  I  error  is  almost  the  same  whereas  the  Type  II  error  increases 
from  the  range  of  13.07  -  19.77  (Table  18)  to  the  range  of  33.57  -  40.37 
(Table  20).  The  Type  II  error  becomes  more  serious  when  the  statistical 
dispersion  of  the  flaw  length  prior  to  inspection  is  smaller. 

Consequently,  it  is  concluded  that  even  with  the  advanced  inspection 
technology  currently  available,  the  Type  II  error  may  be  of  practical 
concern,  depending  on  the  location  of  the  central  portion  of  the  pre¬ 
inspection  flaw  length. 
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Figure  26.  Lognormal  Probability  Density  Function  for 
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TABLE  20 

TYPE  I  AND  TYPE  II  ERRORS  FOR  DISKS  USING  P&W  POD  CURVE  AND 
MEDIAN  FLAW  LENGTH  2.06  mils  (u=0.3133);  LOGNORMAL  DISTRI¬ 
BUTION  FOR  PRE-INSPECTION  FLAW  LENGTH  WITH  50%  DISPERSION, 

+  TOTALLY  CORRELATED  FLAW  LENGTH,  ++  TOTALLY  INDEPENDENT 
FLAW  LENGTH. 


i 

NO.  OF 

HOLES 

^G 

Pl 

^I 

★ 

^I 

1 

92.0% 

3.83x10"^ 

5.20% 

5.65% 

4 

lO'*’ 

92.0% 

2.87x10"^ 

28.60% 

31 .10% 

10^^ 

43.2% 

1.27x10'^ 

19.00% 

44.00% 

1 

97.0% 

8.77x10"^ 

7.36% 

7.59% 

5 

lO'^ 

97.0% 

1.91x10'^ 

33.50% 

34.50% 

lO'^'^ 

73.7% 

3.42x10"^ 

40.30% 

54.70% 

1 

98.8% 

1.89x10"^ 

8.51% 

8.61% 

6 

10^ 

98.8% 

4.00x10"^ 

35 . 30% 

35 . 70% 

lo"" 

88.8% 

7.62x10’^ 

52.70% 

59.30% 

1 

99.5% 

3.74x10"^ 

9.06% 

9.11% 

7 

10^ 

99.5% 

2. 64x10"^^ 

36.00% 

36.20% 

95.3% 

1.52x10’^ 

58.60% 

61.50% 

1 

99.8% 

6.86x10'^ 

9.30% 

9.32% 

8 

10^ 

99.8% 

36.30% 

36.40% 

10^^ 

98.0% 

■witmwfl 

61.20% 

62.40% 

1 

99.9% 

1.16x10'^ 

9.41% 

9.42% 

9 

10^ 

99.9% 

3.56x10"^ 

36.40% 

36.40% 

10^^ 

99.1% 

4.73x10’^ 

62.20% 

62.80% 
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SECTION  VII 

CONCLUSIONS  AND  DISCUSSION 

An  exploratory  study  has  been  made  of  the  possible  application  of 
multiple  inspection  procedures  and  its  potential  pay-offs.  The  solutions 
are  obtained  for  both  independent  NDE  systems  and  dependent  NDE  systems. 
The  resulting  POD  curve  under  multiple  inspections  is  derived  quantita¬ 
tively  from  the  POD  curves  of  individual  NDE  systems  when  the  union  rule, 
the  intersection  rule  and  the  combination  thereof  are  used.  Thus, 
researchers  and  practicing  engineers  can  play  with  different  combinations 
of  NDE  systems  and  procedures  to  arrive  at  an  optimum  strategy  for  their 
particular  purpose. 

Numerical  examples  are  given  using  available  POD  curves  from  certain 
NDE  systems  to  illustrate  the  basic  idea  and  the  application  of  multiple 
inspection  strategy.  It  is  shown  that  multiple  inspections  using  the 
union  rule  alone,  in  general,  reduce  Type  I  error  but  increase  Type  II 
error,  and  the  effect  is  reversed  if  the  intersection  rule  is  employed 
alone.  It  is  further  shown  that  both  Type  I  and  Type  II  errors  can  be 
simultaneously  reduced  significantly  by  the  combined  use  of  union  and 
intersection  rules.  However,  caution  should  be  exercised  in  selecting 
the  third  NDE  system  (or  POD  curve)  to  minimize  possible  adverse  effects. 
The  sequence  of  inspections  to  minimize  the  inspection  cost  has  also 
been  discussed. 

The  multiple  inspection  procedures  presented  include  multiple  inspec¬ 
tions  using  the  same  NDE  system  or  the  identical  redundant  NDE  system  as 
a  special  case,  in  which  the  second  or  subsequent  inspections  are  per¬ 
formed  under  different  inspection  environments.  In  particular,  the  same 
NDE  system  (or  the  identical  redundant  NDE  system)  may  be  used  for  the 
second  inspection  of  the  rejected  holes  of  gas  turbine  engine  disks, 
with  those  holes  being  polished  or  replicated. 

An  exploratory  study  has  been  made  for  the  inspection  reliability  of 
engine  disks  containing  many  holes.  Because  of  mathematical  simplicity, 
two  extreme  cases  for  the  correlation  between  the  flaw  length  in  each 
hole  of  the  same  disk  have  been  studied  in  Sections  IV-V:  (1)  Totally 
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correlated  flaw  length  and  (2)  Totally  independent  flaw  length.  For  the 
case  of  totally  correlated  flaw  length,  the  Type  I  error  is  insignificant 
whereas  the  Type  II  error  may  be  serious  and  therefore  of  practical 
concern,  if  the  bandwidth  of  the  POD  curve  of  the  NDE  system  is  not  narrow 
enough.  For  the  case  of  totally  independent  flaw  length,  not  only  the 
Type  I  error  warrants  special  attention,  but  also  the  Type  II  error  is 
worse  than  that  for  the  case  of  totally  correlated  flaw  length. 

In  reality,  however,  the  exact  solutions  for  Type  I  and  Type  II  errors 
lie  between  these  two  extreme  cases.  Nevertheless,  from  the  NDE  stand¬ 
point,  the  present  results  for  these  two  extreme  cases  provide  important 
information  for  a  greater  insight  into  the  influence  of  the  NDE  inspection 
system  on  RFC/NDE  results.  One  of  the  important  and  significant  conclusions 
obtained  from  the  present  study  is  that  for  the  RFC  system,  the  bandwidth 
of  the  mean  POD  curve  of  the  NDE  system  is  the  most  important  controlling 
factor  rather  than  a  single  value  of  the  flaw  length  associated  with  a 
90%  detection  probability  and  95%  confidence  level. 

The  bandwidth  of  the  POD  curve,  fortunately,  can  be  reduced  to  a 
desirable  level  by  multiple  inspection  procedures  usually  using  two 
inspections  with  the  intersection  rule;  thus  significantly  reducing  the 
Type  II  error.  As  a  result,  for  the  RFC/NDE  system  in  which  both  Type  I 
and  Type  II  errors  are  of  practical  importance,  redundant  NDE  systems 
for  performing  a  second  inspection  of  the  rejected  holes  by  polishment 
or  replication  may  be  worthwhile  to  consider. 

The  present  investigation  allows  for  a  quantitative  description  of 
the  inspection  reliability  of  engine  disks  in  terms  of  Type  I  and  Type 
II  errors.  The  mathematical  solutions  obtained  herein  provide  tools  for 
NDE  engineers  to  manipulate  various  POD  curves  and  the  inspection  limits 
®NDE  oi^der  to  achieve  the  most  beneficial  RFC  inspection  procedures 
for  engine  disks.  Although  the  range  of  the  crack  length  and  POD  curves 
used  in  the  numerical  examples  of  Sections  IV  and  V  are  not  realistic 
for  gas  turbine  engine  disks,  the  trends  and  conclusions  obtained  will 
not  be  altered,  as  explained  in  those  sections. 
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For  the  damage  tolerant  analysis  in  which  the  structural  safety  is 
of  primary  concern,  the  inspection  limit  aj^Q^  has  been  specified  to  be 
the  crack  length  corresponding  to  a  90%  detection  probability  and  95% 
confidence  level.  As  a  result,  many  of  the  POD  curves  appearing  in  the 
literature  are  associated  with  a  95%  confidence  level.  It  should  be 
emphasized,  however,  that  in  the  Retirement-For-Cause  analysis,  the  life 
cycle  cost  (LCC)  is  the  objective  function  to  be  minimized,  and  hence 
both  Type  I  and  Type  II  errors  are  of  practical  importance.  Consequently, 
the  POD  curve  used  in  the  RFC  analysis  should  be  the  mean  POD  curve, 
i.e.,  the  POD  curve  associated  with  a  50%  confidence  level.  If  a  POD 
curve  with  a  95%  confidence  level  is  used,  the  estimation  of  the  Type  I 
error  will  be  too  conservative,  whereas  the  Type  II  error  estimate  will 
be  overly  unconservative.  The  trend  will  be  reversed  if  a  POD  curve 
with  a  low  level  of  confidence  is  employed. 

It  has  been  shown  previously  that  the  distribution  of  the  pre¬ 
inspection  flaw  length  has  a  significant  effect  on  both  Type  I  and  Type 
II  errors.  An  obvious  example  is  that  the  Type  I  error  will  be  zero  if 
all  the  flaw  lengths  prior  to  inspections  are  smaller  than  the  inspection 
limit  In  fact,  the  distribution  of  the  pre-inspection  flaw  length 

is  a  function  of  service  time  {i.e.,  time  dependent),  and  owing  to  crack 
propagation  it  is  shifting  continuously  to  the  large  crack  size  region 
(e.g.,  to  the  right-hand  side  of  Figure  11)  as  the  service  time  increases. 
The  numerical  results  shown  in  this  report  may  represent  the  situation 
for  the  first  few  inspection  maintenances  in  service.  Both  Type  I  and 
Type  II  errors  will  increase  rapidly  as  the  central  portion  of  the 
distribution  of  the  pre-inspection  flaw  length  (i.e.,  a  majority  of  the 
flaw  lengths)  is  shifted  into  the  large  flaw  size  region  close  to  a|^Q£, 
at  a  later  service  time.  Unfortunately,  such  a  situation  is  unavoidable, 
because  the  philosophy  of  the  RFC  system  is  to  exhaust  the  fatigue  life 
of  every  individual  component.  Thus,  during  the  inspection  maintenance 
at  a  later  service  time,  both  Type  I  and  Type  II  errors  will  be  much 
higher  (or  more  serious)  than  those  examples  presented  herein.  As  a 
result,  multiple  inspection  procedures  proposed  may  be  even  more  bene¬ 
ficial  as  the  service  time  increases.  Further  study  is  needed  in  this 
regard. 


AFWAL-TR-82-4111 


While  the  Type  I  error  is  directly  related  to  the  safety  and 
reliability  of  engine  disks  in  service,  the  structural  reliability  can 
be  alleviated  or  safeguarded  by  a  safety  factor  incorporated  in  the 
return  to  service  interval.  Therefore,  there  is  a  remedy  available  for 
the  Type  I  error.  However,  if  the  bandwidth  of  the  POD  curve  is  not 
narrow  enough,  the  only  alternative  to  cope  with  the  Type  II  error,  that 
rejects  good  disks  during  the  inspection  maintenance,  appears  to  be  the 
multiple  inspection  procedures,  in  particular  when  a  fleet  of  disks  is 
in  service  for  a  long  time. 

The  inspection  reliability  of  gas  turbine  engine  disks  using  advanced 
eddy-current  laboratory  inspection  technology  (recently  reported  by  Pratt 
&  Whitney  Aircraft)  has  been  investigated.  Such  an  inspection  system  is 
capable  jf  detecting  a  5  mils  (5x10“^  inches)  flaw  with  a  55%  proba¬ 
bility  and  its  POO  curve  is  assumed  to  have  a  lower  bound  at  1.4  mils 
below  which  no  flaw  can  be  detected.  The  results  for  the  Type  II  error 
may  be  of  practical  concern,  depending  on  the  location  of  the  central 
portion  of  the  distribution  of  the  flaw  length  prior  to  inspection. 

When  the  Type  II  error  is  not  acceptable,  the  multiple  inspection 
procedures  proposed  herein  should  be  considered. 
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APPENDIX  A 

MULTIPLE  INSPECTIONS  WITH  CORRELATED  NDE  SYSTEMS 

Let  R.  be  the  event  that  the  crack  length,  a,  is  rejected  (detected) 
by  the  jth  NDE  system.  Then,  we  have  P0D(a;j)=P[R .]. 

\j 

(1)  Intersection  Rule:  The  probability  of  detecting  the  crack  length, 
a,  by  both  No.  1  and  No.  2  NDE  systems  is  given  by 

P0D(a;lQ2)  =  PCR^aR^]  =  P[R^ ]P[R2l ] 

=  P0D{a;l )P0D(a;2| 1 )  (A-1 ) 

in  which  P0D(a;2|l)  =  conditional  probability  that  the  crack  length,  a, 
is  detected  by  the  No.  2  NDE  system  under  the  condition  that  it  has 
been  detected  by  the  No.  1  NDE  system.  The  experiments  to  be  performed 
to  establish  such  a  conditional  POD  curve  is  described  in  Appendix  B. 

If  both  NDE  systems  are  statistically  independent,  then  P0D(a;2|l)  = 
P0D(a;2)  and  Equation  A-1  reduces  to  Equation  4,  i.e., 

P0D(a;lfi2)  =  P0D{a ;1 )P0D(a  ;2)  (A-2) 

In  a  similar  fashion,  the  resulting  POD  curve  for  m  dependent  NDE 
systems  can  be  derived,  except  that  there  will  involve  more  conditional 
POD  functions. 

(2)  Union  Rule:  The  probability  of  detecting  the  crack  length,  a, 
by  either  the  No.  1  NDE  system  or  the  No.  2  NDE  system,  or  both,  is 
given  by 


P0D(a;lU2)  =  PLR^UR^]  =  PLR^]  +  P(R2]  -  P[RiflR2] 

=  P[Ri]  +  P[R2]  -  P[Ri]P[R2lRi] 


or 

P0D(a;lU2)  =  P0D(a;l)  +  P0D(a;2)  -  P0D{a -,1  )P0D(a  ;2|  1 )  {A-3) 

where  P0D(a;2ll)  has  been  described  before. 
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If  both  NDE  systems  are  statistically  independent,  we  have 
P0D(a;2|l)  =  P0D(a;2)  and  Equation  A-3  reduces  to  Equation  2, 

P0D(a;lU2)  =  P0D(a;l)  +  P0D(a;2)  -  P0D(a  ;1 )P0D(a ;2 )  (A-4) 

Similar  procedures  can  be  applied  to  m  correlated  NDE  systems. 

(3)  Combination  Rule:  The  POD  curve  resulting  from  the  application 
of  three  inspections  with  union-intersection  rule  shown  in  Figure  7  can 
be  expressed  as 

P0D[a;(lU2)fJ3]  =  P[(R^UR2)«R3]  =  P[R3]P[R^  UR^I  R3] 

=  P[R3]{P[Ri|R3]  +  P[R2|R3]  -  P[RiSlR2lR3]} 

=  P[R3]{P[Ri|R3]  +  P[R2|R3]  -  P[Ri  lR3]P[R2|Rif^R3]> 
or 

P0D[a;(lU2)fi3]  =  P0D(a  •,3){P0D(a  ;1  |  3)  +  P0D(a;2l3) 

-  P0D{a;l  |3)P0D[a;2|(la3)]}  (A-5) 

in  which  P0D(a;l|3)  =  conditional  probability  that  the  crack  length,  a, 
is  detected  by  the  No.  1  NDE  system  under  the  condition  that  it  has  been 
detected  by  the  No.  3  NDE  system,  and  P0D[a ;2| (ln3)]  =  conditional 
probability  that  the  crack  length,  a,  is  detected  by  the  No.  2  NDE  system 
under  the  condition  that  it  has  been  detected  by  both  No.  1  and  No.  3 
NDE  systems. 

Usually  the  No.  3  NDE  system  is  required  to  have  a  high  resolution 
capability.  If  the  No.  3  NDE  system  is  statistically  independent  of  No.  2 
and  No.  1  NDE  systems,  then  Equation  A-5  is  simplified  as 

P0D[a;(lU2)n3]  =  P0D(a  ;3){P0D(a ;1 )  +  P0D(a;2) 

-  P0D(a;l )P0D(a;2|l )}  (A-6) 

in  which  P0D(a;2|l)  has  been  explained  before. 
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If  both  No.  1  and  No.  2  NDE  systems  are  also  statistically 
independent,  then  Equation  A-6  reduces  to  Equation  7  as  follows 

P0D[a;(lU2)n3]  =  P0D(a;3){P0D(a;l )  +  P0D(a-,2) 


-  P0D(a;l )P0D(a;2)} 
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APPENDIX  B 

EXPERIMENTAL  PROCEDURES  FOR  ESTABLISHING  CONDITIONAL  POD  CURVE 

The  conditional  probability  of  detection  P0D(a;2|l)  by  the  No.  2  NDE 
system  under  the  condition  that  the  crack  has  been  detected  by  the  No.  1 
NDE  system  can  be  established  by  the  following  procedures.  (1)  N  speci¬ 
mens  with  the  same  crack  length  are  manufactured  and  inspected  by  the 
No.  1  NDE  system.  (2)  After  inspection,  specimens  are  rejected  (or 
detected)  and  N-N-j  specimens  are  accepted  (not  detected).  (3)  Only  the 
N^  specimens  rejected  by  the  No.  1  NDE  system  are  further  inspected  by 
the  No.  2  NDE  system;  with  specimens  being  rejected  and  N^j-N^  speci¬ 
mens  being  accepted  by  the  No.  2  NDE  system.  (4)  The  conditional 
probability,  P0D(a;2|l),  is  the  ratio  of  N^  to  N-j  as  N-h®,  i.e., 
P0D(a;2|l)  =  |||^^(N2/N^  ) .  (5)  Procedures  (1)  through  (4)  are  repeated 

for  different  crack  length.  A  flow  chart  is  shown  schematically  in 
Figure  27. 

The  procedure  described  above  is  straightforward  but  may  be  tedious. 
It  is  expected  that  more  effir  ent  procedures  can  be  established  in  a 
further  research.  -  - 
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Figure  27.  Procedures  for  Establishing  Conditional  POD  Curve. 
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